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INVESTIGATIONS  PERFORMED  ON  THE  ARCTIC  ICE 
DYNAMICS  JOINT  EXPERIMENT,  MARCH  1971 

by 

S.F.  Ackley,  W.D.  Hibler.  Ill,  A.  Kovacs, 

W.F.  Weeks,  A.  Hartwell  and  W.J.  Campbell 


INTRODUCTION 

The  /.retie  Ice  Dynamics  Joint  Experiment  was  formulated  in  the  latter  part  of  the  1960’s 
(Untersteiner  and  Hunkins  1969)*  when,  based  on  previous  measurements  and  models,  it  became 
apparent  mat  a  coordinated,  multidisciplinary,  international  investigation  would  be  useful  for 
determining  the  necessary  parameters  for  a  predictive  model  of  the  motion  and  deformation  of  pack 
ice.  This  predictive  capability  could  then  be  applied  to  such  problems  as  navigation,  climate  pre¬ 
diction,  construction,  and  over-ice  transportation  in  the  Arctic  Basin.  A  revised  scientific  plan 
lias  been  formulated  to  implement  such  a  program  (Maykut  et  al.  1972).  t  The  basis  of  the  plan  is  to 
occupy  a  set  of  stations  (four  in  the  current  plan)  on  the  arctic  pack  for  a  one-year  period  and  per¬ 
form  a  variety  of  meteorological,  oceanographic,  and  ice  dynamics  experiments.  These  experiments 
would  specify  the  driving  forces  and  ice  response  to  the  driving  forces,  thus  allowing  formulation 
of  predictive  physical  and  numerical  models.  This  main  experiment,  currently  scheduled  for  1976, 
was  preceded  in  1971  and  1972  by  a  set  of  pilot  studies  in  which  camps  were  occupied  for  shorter 
periods  of  time.  Efforts  were  made  to  develop  techniques,  test  equipment  and  provide  data  for  the 
formulation  of  the  models  and  to  plan  for  the  main  experiment.  Post-operations  reports  of  these 
studies  appear  in  AIDJEX  Bulletins  No.  7**  and  No.  14. tt  The  pilot  studies  have  proved  to  be 
valuable  experiments  in  t heir  own  right,  yielding  new  data  and  new  interpretations  in  many  different 
areas  of  arctic  sea  ice  investigations,  arctic  meteorology  and  Arctic  Ocean  oceanography. 

This  report,  divided  into  five  parts,  describes  the  work  performed  by  the  USA  CRREL  group 
working  on  AIDJEX  using  data  obtained  during  the  1971  pilot  study.  Parts  1,  2  and  3  have  been 
published  in  the  literature.  Parts  4  and  5  have  not  previously  been  reported.  The  purpose  of  this 
report  is  to  present  in  final  forn  a  single  compilation  of  the  ice  deformation  and  morphology  data 
obtained  during  this  program  i”  1971. 


♦Untersteiner,  N.  and  K.L.  Ilunkins  (1969)  Arctic  Ice  Deformaiion  Joint  Experiment,  Final 
Report.  University  of  Washington,  ONR  Contract  N  00014-67-A  0103-0004,  39  p. 

tMaykut,  G.A.,  A.S.  Thorndike  and  N.  Untersteiner  ( 1373)  AIDJEX  scientific  plan.  Univer¬ 
sity  of  Washington,  Division  of  Marine  Resources,  AIDJEX  Bulletin  No.  15,  67  p. 

‘♦AIDJEX  Bulletin  7  (1971)  1971  Pilot  study  post  operations  narratives.  University  of 
Washington,  Division  of  Marine  Resources. 

++A1DJEX  Bulletin  14  (1972)  1972  Pilot  study  post  operations  narratives.  University  of 
Washington,  Division  of  Marine  Resources. 


PART  I:  MESOSCALE  STRAIN  MEASUREMENTS 
ON  THE  BEAUFORT  SEA  PACK  ICE 

by 

W.D.  Hibler.  HI.  W.F.  Weeks.  S.  Ackley. 

A.  Kovacs  and  W.j.  Campbell 


ABSTRACT 


The  deformation  of  a  strain  triangle  approximately  6  x  8  *  11  km  located  on 
X*  ^  ^  ^  ~  ““  *  -week  period  ZT 

periods  Tie  Ion  T'T'T  (  L5%>  W6re  °bserved  durin*  sl“>«  (6-hour)  tin.e 
n  OR  ‘  ’m-3  ,  -t  cTm  °ne  ay  °r  m0re)  dlver«ence  rate  varied  between  0.04  and 
10-3  hr-l  7  k  :,er"‘  diver«ence  rate*  showed  values  as  high  as  0.29  x 
10  hr  .  The  observed  shearing  motion  indicated  that  the  floes  to  the  east  were 

rela,lve  t0  tlu“  floes  »°  ‘he  west.  This  agrees  with  the  shear  pattern 
that  might  be  expected,  considering  the  location  of  the  station  in  the  Pacific  Gyre 
udies  of  fractur  '  (lead  and  crack)  orientations  in  the  vicinity  of  the  strain  triangle 
indicate  reasonable  correlations  with  the  orientation  of  the  strain  rate  ellipse  ‘  A 
qua  native  relation  is  suggested  between  the  fracture  density  and  the  long-term 

wind  fieTd"  ra,e'  10,18  ^  alS°  °bSerVed  bP,Ween  the  Mergence  of  the 

leld  as  computed  from  the  surface  pressure  field  and  the  ice  divergence. 


INTRODUCTION 

teionesdaand7ber  °, ,lnVeS'1Rat0rs  who  have  a,,fimPted  to  analyze  the  drift  tracks  of  manned  tee^U- 
eset  s  ips  ampbell  1968).  The  first  of  these  investigations  (Nansen  1902)  considered 

equip  ,0'Ce  C  “  "*  ’*""**'  «■  >"  *°n>entum 

boundary  layers.  All  of  these  forces  can  be  considered"^"06 t  ^  PX1Stence  of  aPProPnate 
small  areas  of  interest,  such  as  the  One  nn  .  ®  dered  to  exist  at  every  Poil“  on  the  ice.  For 
Principle,  be  measured  and  expressed  in  terms *7  *  Statl0n  ls  loca,ed-  these  forces  can,  in 
progressing,  field  observations  deal  Even  as  was 

pack.  For  example,  severe  ridging  frequently  oc  *  bl^nif*rant  lateral  transfer  of  stress  through  the 
lateral  transfer  of  stress  has  been ^  ^  ^  winds  were  calm.  This 

forces  currently  included  in  the  momentum  equation  (1  'to  JOT™  '  ™d|  'S  "‘e  le0SI  t,nde,stooli  of  the 
internal  stress  field.  1  V  •  r).  The  „Te  „r  1  T  PreC1Se‘V  «“  <"«««»»  of  'he 
and  the  distribution  of  ice  types  as  well  as  re  h  V  f,lncI10n  of  bo,h  lhe  regional  stress  field 
leads,  in  short,  i,  is  deiormE 'by  the Lltl ’"eSS  e'eme"*s  «  fifties,  hummocks  and 
occur,  the  resultinu  deformation  features  modify  the10  ^  S'raln  hlslorv  of  tlle  ice-  As  few  strains 

surfaces.  This  Cannes  the  vies  ,  Td  f  xerl TthTf ' ”1°'  """  -  “>*«' 

affects  /.  a  w  exerted  on  the  surrounding  pack,  which  in  turn 
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ARCTIC  ICE  DYNAMICS  JOINT  EXPERIMENT 


The  internal  ice  stress  was  first  treated  as  a  simple  frictional  resistance  by  Sverdrup  (1928) 
and  later  as  an  effective  viscosity  by  Ruzin  (1959),  Reed  and  Campbell  (1962)  and  Campbell  (1965). 
The  later  approach  assumes  that  on  a  large  scale  the  ice  pack  can  be  considered  to  act  as  a  thin 
layer  of  a  Newtonian  viscous  fluid.  Using  the  viscous  model  for  /,  Campbell  (1965)  has  been  able 
to  predict  realistic  mean  drift  velocities  as  well  as  the  actual  position  of  the  Pacific  Gyre.  His 
results  also  show  how  sensitive  the  mass  transport  and  ice  flow  fields  are  to  changes  in  /.  Un¬ 
fortunately  his  model  gives  unrealistically  high  convergence  rates  for  the  gyre. 

It  is  possible  to  r  onceive  of  a  number  of  alternative  ways  to  treat  the  interml  ice  stress 
problem  (see,  for  example,  the  papers  in  A1DJEX  Bulletin  No.  2,  1970).  However,  without  an 
experimental  basis  for  evaluating  the  results,  it  would  be  difficult  to  choose  among  the  different 
approaches.  Hence,  one  of  the  most  urgent  needs  is  for  sets  of  good  field  observations  on  the 
actual  deformation  of  the  arctic  ice  pack,  both  under  a  wide  range  of  ice  conditions  and  on  several 
time  and  space  scales.  Because,  in  principle,  it  is  always  possible  to  measure  G,  C,  ra  and  rw 
at  the  deformation  sites,  it  should  at  times  be  possible  to  determine  I  as  a  residual.  The  variations 
in  /  can  then  be  related  to  both  the  regional  stress  and  strain  fields  and  to  differences  in  the  large- 
scale  morphological  characteristics  of  the  pack  in  the  study  areas. 

This  paper  describes  the  results  of  a  pilot  study  of  the  local  deformation  of  the  pack  in  the 
vicinity  of  Camp  200,  the  site  of  the  1971  A1DJEX  pilot  project  in  the  Beaufort  Sea.  The  purpose 
of  this  study  was  to  test  experimental  procedures  and  to  obtain  information  on  both  the  total  strain 
and  the  variation  in  the  strain  rate  that  occurred  during  the  occupancy  of  the  stations. 


PREVIOUS  WORK 

To  the  best  of  our  knowledge,  there  have  been  only  a  few  attempts  to  measure  the  deformation 
of  the  pack  ice  even  though  the  desirability  of  making  such  measurements  has  been  obvious  for 
some  time.  The  reason  for  this  paucity  of  data  is  clear:  most  field  parties  operating  in  the  Arctic 
Ocean  have  been  based  on  only  one  drifting  station  and  have  not  had  the  air  support  necessary  to 
lay  out  and  monitor  strain  arrays  in  areas  of  active  ice  deformation. 

Two  studies  of  local  strains  have  been  made  using  a  theodolite  and  a  base  line  laid  out  ou  an 
ice  islaud.  In  1952  C.ary  followed  the  relative  motion  of  several  prominent  hummocks  located  in 
the  pad;  within  5  km  of  the  edge  of  Fletcher’s  Ice  Island  1-3  (Browne  and  Crary  1958).  The  relative 
locations  of  these  hummocks  did  not  change  significantly  over  a  five-month  period.  In  addition, 
between  1962  and  1965  a  series  of  similar  measurements  was  made  from  ARL1S  11  by  Senior,  Wittmann 
and  Skiles  (1968).  In  the  spring  of  1962,  the  relative  motion  of  four  towers  located  within  1130  m  of 
the  edge  of  the  ice  island  was  monitored  every  two  to  three  days  for  a  one-month  period.  During  this 
time  period  a  lead  9  m  wide  opened  within  the  array  and  then  closed,  forming  a  3-m  pressure  ridge. 

In  the  spring  of  1963,  a  larger  (nine-tower)  array  was  established  and  monitored  at  similar  time 
intervals  for  l'/2  months.  The  maximum  distance  of  a  tower  from  one  end  of  the  base  line  was  just 
under  5  km.  For  the  first  three  weeks,  little  motion  occurred.  Then  a  lead  300  m  wide  rapidly 
opened  between  the  base  line  and  the  towers.  Later,  a  number  of  leads  15  m  wide  developed  in  the 
area  of  the  towers.  Di-ergences  as  large  as  0.039  were  observed,  although  most  values  lay  between 
40.01.  Similar  studies  were  performed  in  the  spring  of  1963  and  in  the  fall  of  1963,  1964  and  1965, 
but  no  significant  relative  motion  was  observed. 

The  T-3  and  ARL1S-I1  studies  share  a  common  difficulty.  Ice  islands  are  not  typical  elements 
of  the  arctic  pack.  In  general,  differences  in  movement  between  ice  islands  and  the  surrounding 
pack  are  to  be  expected.  The  reason  for  this  is  the  difference  in  the  roughness  of  the  upper  and 
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lower  surfaces  of  the  ice  island  as  compared  to  the  ridged  sea  ice,  and  the  larger  Coriolis  effect 
on  the  ice  island  because  of  its  larger  mass.  Sometimes,  however,  the  sea  ice  near  an  ice  island 
appears  to  move  as  a  unit  with  the  island.  For  example,  Senior,  Wittmann  and  Skiles  estimate, 
based.on  observations  from  an  aircraft,  that  during  much  of  the  period  of  their  strain  measurements, 
the  ice  within  8  km  of  ARL1S-1I  was  moving  as  a  unit.  In  short,  unless  independent  verification  is 
available,  strains  measured  in  the  vicinity  of  an  ice  island  axe  not  necessarily  representative  of 
strains  in  the  surrounding  sea  ice. 

Dunbar  and  Wittmann  (1963)  have  also  measured  variations  in  the  areas  of  the  triangle  formed 
by  NP-10,  NP-11  and  ARLIS-ll  and  the  quadrilateral  formed  by  these  three  stations  plus  T-3. 
Changes  of  up  to  20%  were  noted  over  a  15-day  period  (the  measurement  interval),  and  a  maximum 
change  of  63%  relative  to  the  initial  area  was  noted  over  the  total  period  of  measurement  (four 
months).  The  areas  involved  (roughly  140,000  km2  and  280,000  km2,  respectively)  were  latge 
enough  that  differences  in  motion  between  the  ice  islands  T-3  and  ARLIS-ll  and  the  surrounding 
pack  would  not  be  significant.  Other  numerical  measures  of  the  deformation  were  not  calculated. 

The  only  study  of  the  relative  motion  of  a  station  array  located  completely  on  pack  ice  was 
performed  near  the  North  Pole  (approximately  88  N,  147°E)  during  the  spring  of  1961  (Bushuyev 
et  al.  1970).  Four  stations  were  placed  in  a  square  with  sides  of  75  to  100  km.  In  general,  the 
stations  moved  as  a  block,  even  down  to  two  small  counterclockwise  loops  which  occurred  at 
roughly  the  same  time  (within  6  to  8  hours)  in  each  of  the  otherwise  fairly  straight  tracks.  A  study 
of  the  station  velocities  clearly  indicated  that  the  ice  drift  commonly  preceded  the  arrival  of  the 
actual  wind  and  that  the  lead  time  increased  with  the  wind  strength.  Bushuyev  et  al.  (1970)  felt 
that  the  influence  of  the  wind  was  transmitted  through  the  ice  for  at  least  150  to  200  km.  This 
again  emphasizes  the  importance  of  the  internal  ice  stress  term  in  the  momentum  equation.  Al¬ 
though  a  simple  linear  relation  was  shown  between  the  angular  rotation  of  the  array  of  stations  and 
the  vorticity  of  the  wind  field,  no  parameters  relating  to  the  deformation  of  the  ice  were  computed. 

This  survey  of  previous  work  reveals  that  in  a  mesoscale  study  of  ice  deformation  one  might 
expect  strains  of  anywhere  between  zero  and  several  per  cent.  Because  studies  have  not  been 
made  of  the  homogeneity  of  pack  deformation  as  a  function  of  the  array  size,  there  is  little  basis 
for  deciding  the  optimum  size  of  the  strain  array.  Also,  little  is  known  about  how  frequently  the 
deformation  array  should  be  surveyed.  In  the  past  the  most  frequent  observations  were  made  daily. 


SITE  LOCATION 

The  site  of  the  1971  AIDJEX  base  camp  (Camp  2J0)  was  on  the  edge  of  a  large  multiyear  floe 
roughly  16  km  across,  located  at  73°45'N,  130°  15' W  ;Fig.  1).  The  drift  of  the  camp  was  generally 
in  a  south  or  southeasterly  direction  with  the  rate  varying  from  2  to  15  km  per  day. 

Aftei  an  initial  aerial  examination  of  the  pack  in  the  vicinity  of  the  base  camp,  it  was  decided 
to  establish  the  strain  array  in  the  approximately  2-m-thick  first-year  ice  to  the  northeas*.  The 
general  layout  of  the  strain  array  relative  to  Camp  200  and  to  the  major  fractures  and  leads  in  the 
vicinity  is  shown  in  Figure  2.  MRA-3  tellurometers  were  used  as  distance-measuring  instruments. 
Distance  measurements  alo  .g  the  line  a-/3  were  obtained  beginning  on  11  March  with  complete 
triangle  closures  not  being  made  until  12  March.  Several  closures  were  made  thereafter  with  the 
final  closure  being  made  on  23  March. 
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Pigure  1.  Location  of  the  1971  AIDJEX  camp  (position  0).  Also  shown  are  the 
locations  of  the  four  other  positions  where  the  surface  barometric  pressure  was 
obtained  so  that  the  wind  stress  could  he  calculated.  Positions  1,  3  and  4  are 
located  near  permanent  weather  stations.  Pressures  at  positions  1-4  were  taken 
from  Canadian  meteorological  maps.  The  distance  a  (eg  13)  was  about  320  km. 


RESULTS 


Computational  technique 

The  following  procedure  was  used  to  determine  the  strain  rate  tensor  from  the  closure  data. 
Consider  a  medium  with  the  velocity  at  each  point  of  the  medium  given  by  v(x).  Since  we  are  con¬ 
cerned  only  with  the  horizor  tal  motion  of  the  ice  pack,  the  velocity  vector  is  two-dimensional. 
Using  tensor  notation  the  strain  rate  tensor  is  defined  by 


To  measure  the  strain  rate  tensor  we  need  only  to  measure  relative  positions  between  sets  of  points, 
say  P  and  P' ,  as  a  function  of  time.  In  particular,  if  we  choose  a  coordinate  system  such  that  the 
line  PP'  is  oriented  at  an  angle  6  to  the  x-axis,  then  the  one-dimensional  strain  rate  along  PP' , 
denoted  by  r(PP' ),  is  related  to  the  strain  rate  tensor  by  the  equation  (Nye  1957) 

r  j  j  cos^  0  ►  r22  sina  0  r  2<  12  sin#  cos  0 . 


J(PP' ) 


(2) 
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II  Mar.  '71  12  Mar  '71 

Figure  3.  Strain  and  strain  rates  on  the  a-j8  line,  11-12  March  1971. 


This  result  is  easily  obtained  by  transforming  into  a  coordinate  system  with  the  x-axis  parallel 
to  PP' .  By  measuring  linear  strain  rates  along  three  or  more  non-colinear  lines,  a  set  of  equations 
of  the  form  of  eq  2  may  be  solved  for  Since  in  our  case  only  three  linear  strains  were  uaed,  the 
set  of  equations  yielded  unique  values  for  the  strain  rate  tensor  so  that  least  squares  averaging  was 
not  necessary.  (For  a  discussion  of  the  least  squares  solution  of  eq  2,  see  Thorndike  1970.) 

Linear  strains 

As  was  discussed  earlier,  technical  problems  prevented  closing  the  strain  triangle  until 
12  March.  However,  on  the  11th,  an  interesting  set  of  detailed  readings  was  obtained  between  sites 
a  and  )3.  The  results  are  shown  in  Figure  3.  The  total  strain  r  is  given  by 
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ft  \  14  1  19^  16  '  ft  +  iS+i9^  20  »  it  '  U  '  a 
Mar  '71 

Figure  4.  Strain  on  the  a-j9  line,  11-23  March  1971.  The  solid 
portio  is  ol  the  curve  indicate  times  when  the  length  ol  the  line 

was  measured. 


f  = 


A/ 

To 


(3) 


where  i0  is  the  initial  length  of  the  strain  line  (in  this  case  8317  m).  Figure  3  shows  that  after  an 
hour  of  no  variation  in  f,  a  rapid  extension  started  with  expansions  of  as  much  as  0.7  m  min  . 

•  After  4.5  hours,  the  extension  essentially  ceased  (at  1730  hours)  and  tnen  quickly  started  again. 

If  we  examine  the  simultaneous  plot  of  the  strain  rate  i,  we  see  that  it  reached  a  maximum  of  0.116 
day-1  at  1545  hours,  then  decreased  to  roughly  zero  at  1730  hours,  then  rapidly  increased  to  0.142 
day-1  before  finally  decreasing  to  a  roughly  constant  value  of  0.111  day"1.  These  events  suggest 
that  we  saw  the  transfer  of  momentum  through  the  pack  by  jostling  between  floes.  At  the  time  we 
returned  to  the  base  camp,  the  rapid  extension  was  still  continuing.  When  we  reoccupied  the  a-/3 
line  at  0945  hours  the  following  day.  the  value  of  t  was  gradually  approaching  zero,  indicating  a 
near  return  to  the  state  of  strain  that  existed  prior  to  the  extension.  The  maximum  period  of  time 
over  which  this  extension  occurred  was  roughly  22  hours. 


Figure  4  shows  the  strain  along  the  a-f 3  line  for  the  complete  period  of  our  measurements. 
Figures  3  and  4  tell  us  a  great  deal  about  the  necessary  rate  of  data  acquisition  for  mesosoale  strain 
measurements.  Daily  readings  are  clearly  not  adequate:  if  we  had  measured  the  a-j 3  line  at  noon  on 
11  and  12  March,  we  would  have  concluded  that  no  strains  had  occurred  during  this  time  period.  To 
clearly  define  an  event  such  as  this  extension,  readings  should  bt,  taken  at  least  every  15  minutes. 

If  one  wished  to  examine  the  fine  structure  of  the  event,  such  as  the  cusp  in  the  i  curve  at  1730 
hours  on  11  March,  readings  would  be  required  at  least  every  5  minutes. 


Strain  tensor  and  strain  rate  tensor 

The  complete  strain  triangle  was  measured  over  a  10-day  period  beginning  with  an  initial 
closure  at  1100  hours  on  12  March.  Subsequent  closures  were  made  on  15,  21  and  23  March.  On  all 
days  except  23  March,  several  sets  of  closures  were  made  over  a  period  of  several  hours,  allowing 
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th* °flT  h0arS  and  dUyS-  Strain  ratPS  ‘*»P*ed  over 

- ZZZ'TTT  T“  —  ~ 

S==r==SHr 


Table  I.  Net  strains  along  the  legs  of  the  strain  triangle  (units  of  10*3) 

rCi,d,nKM12  1100):  «-A  8324.5  HC-fi,  5579.2  m;  BCw,.  .0.729.0  m 

Date/Time 

(March  1971)  tl.H  nr.H 


12  1100 
12  ir>oo 
if,  iooo 

15/1300 
15  1000 
21/1200 
21/1500 
23/1500 


-0,870 

-0.46 

-1.02 

-1.26 

-2.21 

-2.24 

-2.22 


Jizz  ,6p  d,rec"°"  -  - ««-  « 

lies.  Within  „rls  Hror  h  't  "  ‘  ;  **  “  P"'SS">le  '»  "**"•>  accurate  vorlici. 

.««.  «.  slioi'-teim  strain  tales.  ^  *"*  , 

remained  fixed.  Even  if  this  is  not  ex  u  r'v  r.„»  n  the  a-fi  line  and  true  north 

I*  correct  as  they  arc  indepenZ^ "*  —  *“  ~ 

by  sunm“,,s  ,,,e  m,e™h"  «*«  or  Ihe  strain  icnso,  to  the  ***”* 

of  the  tensor  along  the  pn n c i pal" axe^ (GhTl' (70)" *  t"  ^  ^  a,‘n  ‘  ^  'enS°r  11X6  the  tW0  coniP°nents 
Firmre  5  T.  .  r  1  1  xes  'Ult  "  197°)-  Tlie  results  are  presented  in  this  form  in 

perpendicular  "*£7*!'°*  Klao»»  — >V  as  two 

arc  proportional  to  II, e  mapmudes  of  the  l"W* 

indicates  extension  and  a  bar  denotes  compression  Frnn,  h  ,■  .’  M  ‘‘l‘  outwurd  arrow 
components  of  the  strain  rate  tensor  in  anv  enorrli  \  0lma,I0n  111  tllls  figure,  tlte  shear 

of  the  data,  with  the  exception  of  the  strain  rate  obtc"  v^lbet  ween  p  TdTCh  ^  ,rP"d 

in  approximately  the  east-west  direction.  U  d  lo  M  h’  ls  an  'lxtension 

in  r,  rrr^w  - 

^  ^  — u’5  °"i;;  ;h<> 

'or  later  analysis,  the  dive,Bence!  dlvZncc  , a,,  , ,  7T'’”  *•»  "  "»  » 

presented  in  Figure  6  As  can  be  seen  t.t  1  *  ’  ‘  d  bh(  ai  (V'JXIS  in  the  north  direction)  are 
to  0.08  X  lO-’h,-'  The ,  '*•«  varies  f, on,  0.012  .  10-’  hr"' 

■  U,e  slum  .term  diverponce  rate,  on  tl,e  oilier  Imnd.  Has  values  as  larpe  as 
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f '*“re  5‘  PrinciP‘j/  J«‘s  components  ot  the  strain  rote  tensor  as  a  /unction  ot  time.  The  long-term 
rotes  were  calculated  using  time  intervals  of  two  or  more  days,  whereas  the  other  rates  were  cal¬ 
culated  using  time  intervals  of  Irom  one  to  two  hours.  The  directions  ot  the  bars  indicate  the 
principal  axis  directions  with  their  lengths  being  proportional  to  the  strain  rates 
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Figure  6.  The  divergence  rate  net  divergence  <  and 
net  shear  rxy  as  a  function  of  time.  The  short-term'divcr- 
gence  rites  were  calculated  using  time  intervals  ot  from 
one  to  t ,  o  hours  whereas  the  long-term  rates  were  calcu¬ 
lated  using  time  intervals  of  two  or  more  days. 
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Figure  7.  Net  strain  ellipses  as  a  function  of  time.  The  major 
ellipse  axis  is  in  the  direction  of  maximum  extension  (or  minimum 
compression)  and  the  minor  ellipse  axis  is  in  the  direction  of 
minimum  extension  (or  maximum  compression) 


This  indicates  that  on  the  time  scale  of  hours  large  strains  occur  which  are 


0.293  -  10~3  hr-1 

averaged  out  when  measurements  are  made  over  time  intervals  of  several  days.  The  presence  of 

this  type  of  fluctuation  is  borne  out  by  the  detailed  measurements  along  the  single  strain  line  a-R 
as  shown  in  Figure  4. 

for  a  physical  interpretation  of  the  divergence,  it  is  useful  to  recall  that  when  the  divergence 
rate  is  constant  over  a  region,  the  divergence  rate  equals  the  change  in  area  per  unit  time  divided 
by  he  area,  from  the  plot  of  the  divergence  in  Figure  6.  we  see  that  the  change  in  area  from  12 

the  ke  divergmK  °Ut  ^  ""h  while  from  15  t0  22  March  it  was  about  1%  with 

.  tllA  “SefUl  in  vi&ualizi'V  the  Reformation  in  the  strain  area  during  the  measurement  period 
s  tin  strain  ellipse,  which  is  defined  in  the  principal  axis  system  by 


(1  •  r/  (1  ,  (y)~ 


1  . 


(4) 


Th  s  ellipse  has  a  minor  axis  in  the  direction  of  compression  (or  leas,  extension)  and  a  major  axis 
tn  the  direction  of  extension  for  least  cotnpre-  non).  To  show  the  net  strain  by  using  this  procedure 
wt  have  constructed  net  strain  ellipses  for  15  March  (1300).  21  March  (1200)  and  23  March  0500) 

2En»Ji,"7r  h0urs'  10  «  <*■*>  a.  .tae  zeroon 

‘  ! i's  T  aFe  prespn,ed  in  Fi^ure  7  with  principal  axis  components 

C  y  10°'  1  1S  C  oar  thut  most  of  the  strilin  is  occurring  along  a  line  approximately  in  the 
east-southeast  to  west-northwest  direction. 

Estimation  of  errors 

0  t  °[  the  MRA'3  ,elluron,eters  when  operated  at  only  one  cavity  tune  is  better  than 

0.1  m  for  differential  measurements  of  distance  and  better  than  1  m  for  absolute  measurements.  In 
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view  of  this  type  of  accuracy  an  upper  error  limit  would  be  1.5  m  over  long  time  intervals  (one  or 
more  days  and  0.1  m  over  short  time  intervals  (a  few  hours).  Using  these  values,  the  average 
long-term  linear  strain  measurement  error,  o.  is  0.02%.  and  the  short-term  error  is  0.0013%. 

To  estimate  the  variation  in  the  strain  tensor  components  due  to  such  measurement  errors  we 
will  assume  that  the  errors  are  normally  distributed  with  the  same  variance  along  each  strain  line. 

ith  these  assumplions,  the  maximum  likelihood  (Mathews  and  Walker  1965)  values  (in  the  least 
squares  sense)  for  the  strain  rate  tensor  are  obtained  by  solving  the  set  of  equations  of  ihe  form 
of  eq  2.  The  errors  in  the  strain  rate  tensor  components  within  this  approximation  are  given  by 


Ax  -  a  v/tfLD-M 


mm 

where  Axm  denc'es  the  strain  rate  tensor  errors  with 


(5) 


Ax.  =  A; 


AXj,  -  Afyy. 


Ax„ 


a; 


and 


L(1  =  cos20j,  Lj2  =  sin2^,  Lj3  =  2  sinflj  cos0f  . 

Here  Q.  is  the  angle  between  the  ith  strain  line  and  the  east-west  axis  and  L  is  the  transpose  of  L. 
or  our  particular  array,  using  the  angles  obtained  from  the  first  closure,  eq  5  yields 

A;xx=1.40o,  Aryy  =  0.95  a,  Arxy  =  1.46  a. 

Using  the  values  of  a  mentioned  above  and  eq  5.  we  arrive  at  long-term  net  strain  errors  for  the 
divergence  and  shear  respectively  of  0.047%  and  0.029%.  For  the  strain  rate  errors,  the  short-term 
divergence  and  shear  errors  over  a  three-hour  interval  are  0.01  x  10-3  hr-1  and  0.006  x  10-3  hr-1 
respectively  The  long,  term  divergence  and  shear  errors  over  a  72-hour  interval  aie  respectively 
0.0008  x  10  hr  and  0.005  x  10  3  hr  ^  These  errors  are  presented  as  error  bars  in  Figure  6 

and,  as  can  be  seen,  are  quite  small  compared  to  the  magnitude  of  the  strain  events.  These  errors 
do  not.  of  course,  represent  the  magnitude  of  the  variations  in  the  strain  due  to  lateral  inhomogen- 
eity  in  the  strain  field  (i.e.  different  strains  for  different  triangles).  We  expect  that  such  variations 
would  be  much  larger  than  the  measurement  errors.  A  detailed  study  of  mesoscale  strain  inhomogen- 
eity  is  presently  underway  using  data  collected  in  1972. 


CORRELATION  OF  SYNOPTIC  AERIAL  PHOTOGRAPHY  WITH  MEASURED  STRAINS 

Aerial  photography  of  the  strain  triangle  and  surrounding  area  was  obtained  from  NASA  over¬ 
nights  on  11  and  15  March  and  by  NAVOCEANO  overflights  on  21  and  23  March.  These  dates 
coincided  with  strain  measurements.  Unfortunately  there  was  considerable  variability  in  tho 
photography  because  of  changes  in  tile  weather  conditions,  the  altitudes  flown,  and  the  paths  flown 
Overlays  of  mosaics  of  the  area  from  15  and  23  March  are  shown  in  Figures  2  and  8.  Scale  and  area 
coverage  differences  arise  from  the  altitudes  at  v  hich  the  two  flights  were  made  -  10  600  m  on 
15  March  and  1500  m  on  23  March. 

From  the  photography  or  which  Figure  2  is  based,  net  differences  in  ice  type  could  be  dis¬ 
tinguished  and  are  indicated  on  the  overlay.  These  distinctions  were  based  on  the  surf, cal  appear- 
ance  of  the  different  ice  types.  First-year  ice  is  characterized  by  sharp  angular  ridges,  with 
rela  ively  Hat  ice  between  them.  Multiyear  ice  shows  a  freckled  appearance  caused  by  its  undulating 
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Figure  S.  Overlay  ot  an  infrared  mosaic  ol  the  region  near  Camp  200 
on  23  March  taken  at  an  altitude  ol  1300  m  by  N AVOCEANO.  The 
tract ures,  multiyear  ice  and  thin  annual  ice  were  identified  by  light 
and  dark  tones  on  the  infrared  mosaic. 

■’rm'lt”  topography.  Ice  islands  are  distinguished  by  their  rolling,  washboard-like  topography  and 
wiihin  them  ^  ^  lineS  ,,eCa,,st“  of  ,he  °P''n  wa,er  or  very  thin  ice  present 

Because  Figure  8  was  made  from  imagery  obtained  at  a  lower  level,  more  detail  is  present 
although  the  a*ea  coveted  ,s  not  so  extensive.  This  specific  overlav  was  prepared  from  an  infrared 
mosaic  ,n  uhich  d.Uereiit  thicknesses  of  ice  me  shown  quite  clearly  because  of  differences  in  their 
'  n  lace  temperatures.  I  he  thickest  ice  in  the  mosaic  was  the  multiyear  ice  (2.5  to  5  0  m)  This 
“  <llllk  0,1  th‘‘  III0SiU('  a,,d  is  shc,W11  in  lower  left  of  the  overlay.  Lighter  tones  indi- 

I/ T'1*"  “  2  ""rk)  np  ™st  "»  in  the  area  ol  the  strain 

lunUe  w, thill  this  matrix  of  tirst-year  ice.  two  types  of  thinner  ice  existed.  The  thicker  of  these 

two  mane  tip  the  2000-m-wide  refro/et.  lead  that  crosses  the  BC-n  and  a- ft  lines  and  is  shown  by 
.nPP' '  ‘  *nK  F  itimv  8'  lts  ,,xa(  t  'hickness  was  not  measured,  but  it  is  estimated  to  be  2  1  m 

1  hmnf‘St  ir(!  (0  ,0  40  <:m)  0CCH,red  in  the  fracture  systems  that  were  active  during  the  strain 
measuieinents  on  21  and  2(3  March.  This  ice  appears  as  bright,  sharp  lines  on  the  infrared  mosaic 
(highest  sitriace  temperature)  and  is  shown  by  the  dashed  lines  in  Figure  8. 

From  Figures  2  and  8  it  is  apparent  that  the  active  leads  are  no.  strongly  affected  by  ice 
h  kness  variations,  since  on  both  overlays  leads  cut  through  first-  and  multiyear  ,ce  without 
.  lection  One  can  also  see  in  Figure  8  that  the  active  lead  system  does  not  coincide  with  pre- 
show  T  1,1  T  MeaS’  p0SSlt,ly  because  the  orientations  are  slightly  different.  The  thin  ice  area 

s  0,1 11  March  -  -  p"’si-ahly  ««• ■ r°iai'v«>y zz  r 

Details  ol  the  fracture  structure  were  analyzed  by  determining  the  total  fracture  length  in  anv 
and  1  ir“  haf  h  ,r;K'MirL‘  was  broken  into  straight  line  segments  from  100  to  300  n.  in  length 
an(1  lenpt"  —  ««■«>■  The  total  length  in  each  orientation  (10  mimvlt) 
fracture  l’hc  result  tl,en. t*,lmIate<!  and  converted  to  a  percentage  of  the  total  length  of 
cmen,  'bv  I  ,  Sh°'V“  1,1  9-  1,1  addi,10'b  total  length  was  divided  by  the  area 

.  P  ,  •  .  '“p',y  '»  **"•  *  *«■«“«  dl'ns"y  •m.ilosuns  to  the  ndsc  deiisilv  used  l,v  Mork 

ll.ulwol]  *,  ll.ltlo,  ,19721  iro,„aBo  d0„s„y  used  in  hy(l,olut.y  1,1  Mo,  k- 

tJtzrzz  :r::,,nw  -  t**  *  «■*•  and  «•«««»>. «» «-*««»  ^ 

i  souk  c  arc.  Vie  can  hypothesize,  however,  that  the  net  orientation  of 
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Figure  .9.  Histogram  identifying  the  fracture  orientations  obtained  from 
aerial  imagery  measured  in  a  clockwise  direction  from  the  east-west  line. 
Changes  in  fracture  density  agree  with  convergence  or  divergence  of  the 
pack,  indicated  by  the  long-term  strain  rates.  Net  fracture  vectors  on  the 
15th  and  21st  agree  with  the  direction  expected  from  the  short-term  strain 

rates  (Fig.  10). 


the  fractures  (if  any)  should  be  correlated  with  the  state  of  stress  of  the  ice  pack  as  manifested  by 
the  short-term  strain  rate.  To  make  this  comparison,  a  net  fracture  vector  was  calculated  from  the 
results  in  Figure  9  by  using  Krumbein’s  method  as  described  in  Mock,  Hartwell  and  Hibler  (1972). 
This  method  uses  twice  the  angle  of  each  fracture  orientation  to  convert  the  180°  distribution  to  a 
nonsymmetric  360'  distribution. 

We  calculated 


6 


-  arctan 
2 


£  f  sin  2 0  "I 

1 1  cos  26  J 


(6) 


and 


r  [(E  t  sin  2 0)2  ♦  (i.  {  cos  2 0)2]'4  . 


(7) 
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b.  21  Morch  *71 

Figure  10.  Strain  rate  ellipses  (or  short-term  strain  rates  on  15 
and  21  March  1971. 


Here  0  is  the  estimate  of  the  resultant  vector  direction,  r  is  the  resultant  magnitude  of  this  vector 
in  percent,  and  f  is  the  percentage  of  the  total  lead  length  in  each  orientation  category.  For  a 
directionally  random  population,  r  is  zero.  Referring  to  Curray  (1956),  who  provides  a  chart  for 
evaluating  significant  deviations  from  zero  for  r,  given  the  total  number  of  samples,  we  find  that  in 
all  cases  illustrated  in  Figure  9  we  may  reject  the  hypothesis  that  our  sample  is  directionally 
random  at  the  0.01%  level.  The  directions  of  the  resultant  vectors  are  shown  by  solid  arrows  in 


On  two  days,  15  and  21  March,  both  short-term  strain  data  and  aerial  imagery  were  available, 
or  purposes  of  correlation  it  is  useful  to  represent  the  strain  on  these  days  by  the  strain  rate 
ellipse,  which  is  completely  analogous  to  the  strain  ellipse  presented  in  eq  4  except  that  strain 

in^eT65?  °f  S‘rains  1116  USed'  As  in  the  strain  ellipse,  the  strain  rate  ellipse  has  a  minor  axis 
n  the  direction  of  compression  (or  least  extension)  and  a  major  axis  in  the  direction  of  extension 
(or  least  compression).  The  strain  rate  ellipses  for  15  and  21  March  are  illustrated  in  Figure  10 
n  angular  relationship  between  the  major  axis  of  the  strain  rate  ellipse  and  the  net  lead  direction 
may  hypothesized  as  follows.  Once  leads  are  formed,  we  assume  that  the  deformation  of  the 

tW°  m°TS'  Shear  al0ns  the  leadS  and  extension  (or  compression)  perpendicular 
°o  ,  leads-  Consequently.  tor  a  diverging  strain  the  major  axis  of  the  ellipse  would  vary  from 
90  to  45  relative  to  the  lead  orientation,  with  90°  for  pure  extension  and  45°  for  pure  shear.  Con- 
rse  y,  for  a  converging  strain  the  major  axis  would  vary  from  0°  to  45°  relative  to  the  lead 

nnpT  m"’  0  f<*  *  m  C°mpreSSi0n  311(1  45°  for  pure  shear.  Using  this  model,  we  note  that 

2  March  (dlVerg,nR  rate)  the  fracture  orientation  was  at  9°  (from  true  north)  and  the  major  ellipse 
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a.  11  Marcn  1971.  b.  23  March  1971. 

Figure  11.  Sequential  aerial  photography  showing  shear  deformation  along  the  lead  that  crosses 
the  a-/3  line  approximately  1.5  km  from  f3  and  runs  NNW-SSE  (see  Fig.  8). 


axis  was  at  68°,  with  the  difference  being  between  45°  and  90°  as  expected.  Conversely,  on 
15  March,  when  we  observed  a  converging  strain,  the  lead  orientation  was  at  70°  and  the  major 
ellipse  axis  was  at  95°,  with  the  difference  lying  between  0°  and  45°  as  expected. 

Although  these  agreements  are  qualitative,  they  do  point  to  a  possible  way  to  estimate  the 
direction  (within  45°)  of  the  major  strain  axis  from  observations  of  the  fracture  structure  obtained 
by  aerial  imagery.  This,  of  course,  assumes  that  it  is  possible  to  determine  from  the  imagery 
whether  the  ice  pack  is  converging  or  diverging. 

Another  correlation  migut  be  expected  between  th;  long-term  divergence  rate  and  the  lead 
density.  In  particular,  if  between  two  days  the  lead  density  increased,  we  would  expect  that  a 
diverging  strain  had  occurred.  A  decrease  in  lead  density  would  indicate  the  converse.  To  test 
this  hypothesis,  we  have  calculated  the  changes  in  lead  density  and  compared  them  with  the  long¬ 
term  divergence  rate.  The  results  are  summarized  in  Table  II.  As  can  be  seen,  there  is  general 
agreement,  except  for  the  period  21-23  March.  We  note,  however,  that  the  deformation  during  this 
period  was  dominated  by  a  large  shear  along  a  single  north-south  fracture.  This  fracture  can  be 
seen  in  Figure  11  (see  also  Hartwell  1972,  for  a  detailed  catalogue  of  motions  along  leads. 


Table  D.  Observed  changes  In  lead  density  and  the  lon;-tenn  divergence 
rate  of  the  ice  pack  during  the  same  time  period. 

Change  in  lead  density  Divergence  rale 

Time  period  (km~')  < long-term) 


1 1  to  15  March 
15  to  21  March 
21  to  23  March 


1.025  to  0.287  (neg) 
0.287  to  1.259  (pos) 
1.259  to  0.404  (neg) 


Negative 

Positive 

Positive 
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including  this  particular  shear).  This  may  he  a  case  in  which  the  strain  triangle  was  not  averaging 
deformation  over  a  large  enough  area  to  yield  correct  results.  It  is  quite  easy  to  construct  an  ex¬ 
ample  of  a  triangle  with  only  shear  occurring  along  a  single  fracture  which  cuts  through  the  tri¬ 
angle.  Such  an  example  would  yield  a  positive  divergence  where  no  divergence  hud  occurred.  The 
wind  stress  data  discussed  in  the  next  section  also  indicate  that  a  convergence  would  have  bpen 
expected  between  21  and  23  March,  so  that  a  measurement  anomaly  may  well  have  occurred. 

Aerial  imagery  is  useful  for  documenting  the  particular  type  of  shearing  motion.  The  measured 
net  shear  is  indicated  in  Figure  6.  From  the  measurements  it  is  not  possible  to  determine  the  exact 
type  of  slippage  (it  could  be  east-west  or  north-south).  The  en-echelon  fractures  to  the  north  of 
Site  £  (Fig..  8)  suggest  ihat  eastern  floes  are  moving  south  relative  to  the  western  floes.  This 
conclusion  is  verified  by  a  detailed  examination  of  the  sequential  photography  (Hartwell  1972). 

Such  a  net  shear  agrees  with  the  general  motion  of  the  Pacific  Gyre  vat  the  measurements  were 
taken  on  the  edge  of  the  gyre  wuh  the  outer  part  moving  more  rapidly  tnan  the  inner  part,  causing 
a  north-south  slippage  of  the  type  observed. 

In  summary,  the  experience  of  the  1971  pilot  study  has  made  it  clear  that  sequential  aerial 
imagery  is  a  great  help  in  interpreting  strain  observations  in  the  pack  ice.  However,  because 
significant  large-scale  strain  events  occur  during  time  periods  of  the  order  of  hours  (e.g.  the  events 
of  11  to  12  March),  it  is  quite  easy  to  miss  these  events  on  synoptic  aerial  imagery  taken  at  daily 
or  greater  intervals.  It  is  also  clear  that,  at  present,  imagery  cannot  provide  the  details  needed  for 
short-term  correlations  with  meteorological  factors,  nor  with  the  act  al  measurements  of  divergence, 
net  strains  and,  most  important,  strain  rates  that  are  required  for  modeling  purposes.  At  the  present 
time,  ground  stations  must  still  be  relied  upon  for  such  critical  ice  dynamics  information. 


CORRELATION  OF  ESTIMATED  WIND  STRESS  AND  STRAIN 

To  estimate  the  wind  stress  field,  the  wind  velocity  field  was  calculated  from  barometric 
pressure  readings  by  using  the  expression 


u 


1  pi  dy  2  dx 


(8) 


v 


+  K 


1_  dP 
1  p(  dx 


fiP 

dy 


(9) 


where  u  and  v  are  the  x  and  y  components  of  the  wind  velocity,  P  is  the  surface  barometric  pressure, 
P  the  air  density,  /  the  Coriolis  parameter,  and  Kj  and  Kg  are  positive  constants.  The  first  terms 
in  eq  8  and  9  represent  the  geostrophic  wind  with  a  modification  expressed  by  the  constant  due 
to  the  effect  of  the  frictional  drag  along  the  earth’s  surface.  The  frictional  drag  also  causes  the 
surface  wind  to  have  a  component  moving  across  the  pressure  isobars  from  high  to  low  pressure 
(Pet'erssen  1969,  p.  158).  This  component  is  represented  by  the  second  term  in  eq  8  and  9.  The 
exact  values  of  Kx  and  K2  may  vary  depending  on  the  surface  roughness.  Here  we  are  only  con¬ 
cerned  with  qualitative  correlations  so  the  exact  magnitudes  are  not  critical.  The  divergence  of  the 
wind  velocity  field  by  the  above  equations  is  given  by 


V  ■  v  -  Kz  y2P. 


(10) 
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Because  of  the  approximate  form  of  the  barometric  pressure  data  available  to  us,  it  was  impossible 
to  calculate  meaningful  shear  components  of  the  wind  stress  field;  we  therefore  examined  only 
correlations  between  the  divergence  of  the  ice  and  the  divergence  of  the  wind  velocity  field. 

Assuming  the  wind  stress  on  the  ice  is  proportional  to  the  wind  velocity,  a  compressional 
stress  will  be  introduced  into  the  ice  which  is  proportional  to  V  •  v.  This  may  easily  be  see  n  by 
considering  adjacent  infinitesimal  sections  of  the  ice  and  determining  the  force  one  section  exerts 
on  another  due  to  the  wind  velocity  gradient.  With  a  constant  velocity,  for  example,  there  would  be 
no  compressive  stress  because  the  external  shearing  force  would  be  the  same  everywhere  (neglect¬ 
ing  variations  in  surface  roughness). 

To  the  extent  that  other  external  stresses,  such  as  water  stress,  may  be  neglected,  we  would 
expect  the  convergence  rate  to  scale  with  the  compressional  wind  stress: 


ijj  =  -  K(V°~P) 


(11) 


where  K  is  some  parameter  depending  on  the  compactness.  Intuitively  this  result  might  be  expected 
to  apply  to  very  compact  ice  (high  internal  ice  stress)  which  cannot  move  rapidly,  so  that  Coriolis 
and  water  drag  forces  can  be  neglected. 

The  other  extreme  approximation  that  could  be  made  would  be  to  neglect  the  internal  ice  stress 
completely  and  consider,  for  example,  unconstrained  wind-driven  ice  on  the  open  ocean.  In  this 
case  certain  theories  yield  an  ice  divergence  rate  proportional  to  V2P  (with  a  sign  opposite  to  that 
in  eq  11).*  Consequently,  for  purposes  of  empirical  correlation  V2P  is  clearly  a  useful  parameter 
to  compare  with  the  divergence  rate,  with  the  sign  of  the  correlation  (if  any)  giving  some  indication 
of  the  internal  ice  stress  magnitude. 

To  calculate  V2  P,  we  utilized  the  surface  barometeric  pressure  at  the  five  locations  indicated 
in  Figure  1.  With  reference  to  this  figure,  pressure  PQ  (at  Camp  200)  was  provided  by  Hans 
Prelkkinen  of  the  Canadian  Polar  Continental  Shelf  Project.  These  pressures  were  taken  two  to 
three  times  daily  and  were  therefore  extrapolated  to  six-hour  intervals  for  comparison  with  other 
data.  The  position  of  pressure  P4  '  as  taken  to  coincide  with  the  location  of  Sachs  Harbor  on  Banks 
Island,  where  surface  pressures  were  available  at  six-hour  intervals  through  the  Canadian  Meteoro¬ 
logical  Service.  Other  pressures  were  taken  from  surface  isobaric  charts  supplied  by  the  Canadian 
Meteorological  Service,  with  stations  1  and  3  nearly  coinciding,  respectively,  with  the  Mould  Bay 
weather  station  on  Prince  Patrick  Island  and  the  Barter  Island  station  on  the  Alaskan  coast.  Data 
were  available  at  six-hour  intervals  and  a  time  series  of  pressure  data  was  constructed  beginning 
at  1800,  7  March.  All  data  were  converted  to  GMT  minus  six  hours,  which  was  approximately  the 
local  time  at  Camp  200.  The  distance  from  P  to  the  other  pressure  locations  was  about  320  km. 

From  these  five  pressure  values,  the  Laplacian  of  the  pressure  field  was  calculated  by  using  a 
finite  difference  grid,  with  a  the  distance  between  the  grid  points: 


y^P  = 


Pt  +  P2  +  P3  ♦  P4  -  4 P0 


(13) 


Since  here  we  are  only  concerned  with  examining  the  proportionality,  we  took  a/2  =  1. 

The  results  are  presented  in  Figure  12  where  we  have  plotted  -\^P  versus  time,  together  with 
the  divergence  and  divergence  rate  as  functions  of  time.  To  smooth  the  pressure  data,  we  removed 
the  high  frequencies  with  wavelengths  shorter  than  24  hours  by  passing  a  digital,  low  pass,  unity 
gain  filter  over  the  -V  P  curve.  The  filter  has  a  pass  band  extending  from  0  to  (1/96)  hr-1  and  a 
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More-  1971 


ilime  12.  The  Laplacian  ot  the  surface  pressure  field,  the  diver¬ 
gence  rate  and  net  divergence  ns  functions  ot  time.  The  unit  length 
referred  to  is  about  ISO  km.  The  dots  on  the  divergence  curve  repre¬ 
sent  actual  measurements,  with  the  smooth  -.urve  simply  connecting 
the  measurements.  The  estimated  divergence  between  11  and  12 
March  was  based  upon  visual  ground  observations,  atrial  photography 
and  the  linear  strain  along  the  a-(i  line. 


stop  band  beginning  at  (1/24)  hr*' .  The  filter  *as  designed  according  to. the  procedure  discussed 
V  i  tier  (19/1)  and  had  less  than  0.6"£  side  lobe  error.  Also  note  that  in  Figure  12  an  extrapolated 
uvergence  curve  is  given  for  the  time  period  between  11  and  12  March.  This  extrapolation  is  based 
on  observations  of  the  lead  stiueture  by  using  aerial  photographs  as  well  as  ground  observations 
ant.  is  supported  by  the  strain  data  ta’  en  along  the  a-fl  line  during  the  same  time  period. 

As  can  be  seen,  a  reasonable  correlation  exists  between  the  wind  velocity  field  divergence 
and  the  ice  divergence  rate.  The  ice  between  10  and  12  March  indicated  a  divergence  followed  bv 
a  t  onvergence  which  agrees  with  the  estimated  wind  stress  which  is  positive  on  the  11th  and 
negative  on  the  12th.  The  trends  of  the  later  data  after  15  March  are  also  in  agreement,  i  e  the 
wind  stress  divergence  generally  increases  and  the  ice  divergence  rate  is  positive. 

In  general,  the  results  presented  in  Figure  12  suggest  a  correlation  between  the  negative 
^aplacian  of  the  pressure  field  and  the  ice  divergence  rate,  although  the  data  collected  are  not 
sufficient  to  permit  us  to  make  any  definitive  statements  about  such  a  correlation.  This  is 

w!ndCstaresstrUr.d"eif0  °f  th°  ice  m0,10n  “lld  ,he  ne^lectin«  of  other  stresses  besides 

mUy  noted.  ^  dUta  are  nm,ed  t0  find  0111  whether  such  correlations  are  gen- 
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CONCLUSIONS 

A  number  of  conclusions  can  be  drawn  from  our  study  that  have  considerable  bearing  on  the 
A1DJEX  Project.  These  are  as  follows. 

1.  The  analysis  of  mesoscale  strains  in  the  ice  pack  lias  proved  to  be  of  considerable 
interest.  However,  for  such  a  study  to  be  of  maximum  usefulness  in  the  analysis  of  ice  deforma¬ 
tion,  the  strains  should  be  determined  at  some  fixed  time  interval.  The  measurements  along  strain 
line  a-ft  on  11-12  March  suggest  that  this  time  interval  should  be  no  greater  than  one  hour.  When 
significant  ice  deformation  is  occurring,  an  interval  of  15  minutes  would  be  even  more  satisfactory. 
>Ve  do  not  feel  that  tellurometers  are  a  suitable  instrument  for  such  measurements  unless  they  are 
automated  so  that  only  one  end  of  a  strain  line  has  to  be  occupied. 

Therefore,  during  the  1972  AIDJEX  field  program  we  used  both  continuous-wave  and  pulsed- 
laser  range  finders  to  determine  the  distances  between  the  main  camp  and  a  series  of  remote  towers. 

2.  It  is  important  that  strain  measurements  be  made  on  several  scales  so  that  small  scale  floe- 
floe  interactions  as  well  as  larger  scale  “continuous"  deformation  of  the  pack  can  be  analyzed  and 
related.  Only  when  such  studies  are  completed  can  sound  statements  be  made  about  optimum  sizes 
for  strain  arrays.  An  examination  of  the  aerial  imagery  taken  during  the  overflights  (sec,  for  ex¬ 
ample,  Figures  2  and  8)  as  well  as  visual  surface  observations  definitely  suggest  that  in  most  cases 
the  ice  divergence  as  measured  witlnr.  Hie  mesoscale  array  does  indeed  correlate  with  the  ice 
divergence  as  observed  over  a  larger  area.  Perhaps  such  mesoscale  strains  will  serve  as  an  ade¬ 
quate  index  of  the  macroscale  deformation  of  the  pack. 

3.  The  correlations  between  fracture  orientation  and  the  orientation  of  the  strain-rate  ellipse, 
as  well  as  between  the  lead  density  and  the  divergence  rate,  are  encouraging.  They  suggest  that 
aerial  imagery  may  eventually  produce  useful  quantitative  information  on  the  deformation  of  the  ice 
pack.  However,  at  the  present  time  more  combined  ground  truth  and  remote  sensing  data  will  be 
required  both  to  increase  our  confidence  in  such  correlations  and  to  make  them  more  quantitative. 

4.  In  designing  the  present  experiment  it  was  always  hoped  that  the  observed  correlation  be¬ 
tween  the  divergence  of  the  wind  field  and  the  behavior  of  the  ice  would  prove  to  be  simple.  The 
results  are  certainly  encouraging  in  this  regard.  They  suggest  that  detailed  mesoscale  studies  of 
the  relations  between  the  wind  flow  characteristics  and  the  two-dimensional  deformation  of  the  ice 
pack  will  prove  to  be  most  profitable. 
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PART  II:  STRUCTURE  OF  A  MULTIYEAR  PRESSURE  RIDGE 

by 

A.  Kovacs,  W.F.  Weeks,  S.  Ackley  and  W.D.  Hibler,  Ill 

ABSTRACT 


Three  transverse  profiles  across  a  large  pressure  ridge  located  in  the  Beaufort 
Sea  are  presented.  The  ridge  sail  extended  4  m  above  sea  level  and  the  ridge  keel 
13  m  below.  The  cross  sections  of  the  ridge  keel  can  be  described  as  roughly  semi¬ 
circular,  suggesting  that  form  drag  coefficients  for  flow  transverse  to  the  long  axes 
of  multiyear  ridges  may  be  as  high  as  0.8.  Examination  of  several  salinity,  tempera¬ 
ture  and  brine  volume  profiles  shows  that  much  of  the  ice  in  the  ridge  has  a  very 
low  salinity  and  is  quite  strong.  All  the  inter-block  voids  that  initially  existed  in 
the  ridge  at  the  time  of  its  formation  have  been  completely  filled  with  ice.  These 
observations,  coupled  with  icebreaking  experience,  indicate  that  multiyear  ridges 
are,  indeed,  significant  obstacles  to  even  the  largest  icebreaking  ship  and  should 
be  avoided  if  possible.  A  very  large  first-year  ridge  with  a  sail  height  of  12.8  m 
is  also  described.  This  is  the  largest  free  Boating  ridge  yc*  recorded. 


INTRODUCTION 

A  major  problem  in  the  development  of  the  Arctic  is  the  lack  of  a  cheap,  reliable,  large  volume 
transportation  system  capable  of  year-round  operation.  In  many  arctic  areas  such  as  the  Canadian 
Archipelago,  a  sea  route  seems  to  be  the  only  possible  solution.  The  principal  obstacle  to  the 
development  of  such  a.  rouie  is,  of  course,  the  highly  variable  and  at  times  extremely  formidable 
sea  ice  cover.  This  ice  cover  has  restricted  the  operation  of  surface  ships  to  a  few  months  each 
summer.  Even  then  there  are  areas  where  icebreaker  assistance  is  required.  However,  the  recent 
test  cruises  of  the  S.S.  Manhattan  clearly  showed  that  “super”  cargo  vessels  can  be  designed  and 
built  within  the  framework  of  current  marine  technology  that  are  capable  of  winter  navigation  in 
arctic  waters.  To  optimize  any  such  design  one  must  fully  understand  the  nature  of  the  more 
significant  obstacles  that  can  be  expected  in  the  ice  cover.  Similar  information  is  also  needed  for 
the  safe  design  of  offshore  structures  and  harbor  facilities.  Most  such  structures  are  fixed,  and 
cannot  maneuver  around  the  more  obvious  obstacles  but  must  take  the  ice  as  it  comes. 

With  the  exception  of  ice  islands,  which  are  rare  and  can  usually  be  easily  recognized  by  radar, 
the  most  formidable  obstacles  in  the  ice  pack  are  the  pressure  ridges  and  hummocks,  i.e.  those 
linear  to  irregular  accumulations  of  ice  caused  by  the  compressive  and  shear  interactions  between 
ice  floes;  such  features  may  become  very  large  indeed.  The  highest  free-floating  ridge  sail  that  we 
are  aware  of  was  12.8  m  above  sea  level  whereas  the  deepest  keel  recorded  in  submarine  sonar  data 
extended  some  47  m  below  sea  level  (Lyon,  personal  communication).  Impressive  as  these  figures 
are,  ridges  of  this  size  fortunately  appear  to  be  quite  rare.  Indeed,  ridge  sails  in  excess  of  5  m 
are  rare. 

Operational  experience  has  shown  that  the  bulk  properties  of  geometrically  similar  ridges  may 
vary  greatly.  For  instance  the  S.S.  Manhattan  was  able  to  progress  with  relative  ease  through  most 
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first-year  ridges  but  encountered  considerable  HiffiPi.ir.,  -,i 

Actually  it  might  be  expected  that  the  properties  of  first-year  ^re"tly .®lmilar  nm,tiyear  ridges, 
variable.  This  would  primarily  depend  upon  the  decree  m  !  "dKeS  W°Uld  be 

that  make  up  the  ridge  The  best  bonding  should  d  i  ^  developod  betwe«»  the  ice  blocks 

r^dge  is  very  cold.  Many  of  the  ^  T  °*  ™  *"*  1,1,0  ,h*' 

(Weeks  et  al.  1971).  In  a  multiyear  ridge  on  the  other  I  '  a  7  ^  de,a1'  W6re  P00rly  l,0nded 

ean  drain  downward  into  the  core  ^'hc  r  d7e  meUWater  Pr0dW*d  the  sum- 

blocks  rose.be,  i»,„  a  solid,  resist,  11  '  fnu.  “T.*‘  «*«*  "*«“>•  «—«*  .bo  ice 
usually  required  for  lire  voids  between  the  ice  blocks'll!'!  *  A'  leaSI  ,wo  ni°"  SIM|sons  are 

-  «  . . ".year  ridye  was  «»..lj  ™M*  Tta 

by  Russian  iiivestiyators  (Burke  1910.  Zubov  1M5).  *  *  was  'ccoKinzcci  quite  early 

■be  —  *  ,he  ~  ^  obstacles  in 

[lie  state  of  the  ice  fortnmy  them.  The  only  oubli  '  atVJ,10ll.'s  01  eltller  >berr  overall  structure  or 
we  are  aware  is  by  Spichkin  (Gakkel*  1959)  wbo  d’ ‘ '  ,nvest,Rat,on  of  a  multiyear  ridge  of  which 
ridge  where  it  was  intersected  by  a  lead  C  rev  777  '7  CrOSS*sectio"a»  profile  of  such  a 
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Shelf  Project,  the  base  of  tl*  1971 1  ctk  Ice  Dv  ^  t  °f  the  Canadian  P°l«r  Continental 

The  observations  include  the  i  “  ExPe,ln*nt  (A1DJEX>  p^t  study, 

us.ng  standard  surveying  techniques  the  prof  ling  77  T*l  °\ thR  Und  ^  S"0W  thickness 
by  Kovacs  (1971).  and  the  sampling  ^  “  S°"ar  “*h"*ue  utilized 
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profiles 
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Figure  1.  Aerial  view  ol  the  multiyear  ridge  studied  (center  of  photograph )  and  the  surrounding 
terrain.  The  photograph  was  taken  after  the  study  was  completed  and  a  storm  hud  swept  the  ice 
surface  clean  of  a  major  portion  of  its  snow  cover.  First-year  ice  can  be  seen  in  the  upper  left- 
hand  corner.  The  position  of  each  profile  line  is  marked. 


between  the  ice  thickness  obtained  by  drilling  and  by  sonar.  We  have  no  explanation  for  this  except 
to  observe  that  the  drilled  thickness  is  a  point  measurement  whereas  sonar  “samples”  an  apprecia¬ 
ble  area.  In  a  general  way  the  transverse  profiles  of  the  ridge  keel  can  be  described  as  roughly 
semicircular  to  semielliptical.  Drag  coefficients  for  such  form  roughness  elements  are  discussed  in 
Hoerner  (1965)  and  range  from  0.8  for  a  semicircular  shape  to  0.5  for  a  semiclliptieal  shape.  Be¬ 
cause  of  the  “streamlining”  effects  of  the  summer  melt  season  and  the  snowdrifts,  the  drag  coeffi¬ 
cient  for  the  ridge  sail  would  be  significantly  lower. 
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Figure  2.  Transverse  and  longitudinal  cross  sections  o(  the 
multiyear  ridge. 


In  most  of  the  profiles,  the  upper  ice  surface  at  the  edges  of  the  ridge  sails  shows  a  downward 

?o?nCl1011'  Th'S  phenomenon  has  teen  observed  in  most  fust-year  ridges  studied  (Weeks  and  Kovacs 
19/0,  Kovacs  1971).  In  the  case  of  first-year  ridges  this  deflection  is  clearly  caused  by  load  trans- 
Ter  from  the  sail  to  the  surrounding  plate  ice  as  a  result  of  an  isostatic  imbalance  produced  during 
lie  initial  ice  pile-up.  In  the  current  ridge  the  deflection  is  believed  to  be  partly  related  to  snow 
loading  and  partly  to  keel  ablation.  The  latter  would  cause  the  ridge  to  subside  m  the  process  of 
seeking  isostatic  equil.bnum.  It  is  expected  that  during  the  next  melt  season  these  depressions 
will  fill  with  meltwater  fromthe  ablating  ridge  sail.  When  these  ponds  refreexe  the  tollowing  fall 
these  deflected  areas  will  probably  no  longer  be  discernible. 

The  slope  of  the  ridge  sail  was  determined  by  measuring  the  angle  of  the  profile  lines  shown 
in  the  cross-sections  (Fig.  2).  The  slope  varied  between  12°  and  28°  and  averaged  19°.  This  is 
o  less  than  the  average  surface  slope  determined  by  Kovacs  (1971)  for  first-year  pressure  ridges. 
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Figure  3.  Fence  diagram  showing  the  relative  positions  of  cross  sections  A,  B  and  C. 

Tins  difference,  which  may  well  be  general,  is  undoubtedly  due  to  the  modification  of  the  shape  of 
the  ridge  sail  by  ablation.  No  attempt  was  made  to  determine  the  slope  of  the  keel  which  has  a 
bowl-like  shape. 

The  fact  that  pressure  ridges  act  as  snow  fences  is  again  clearly  illustrated  in  the  cress 
sections.  The  greatest  snow  accumulation  was  in  excess  of  1  m  and  occurred  along  the  toe  of  the 
sail. 


INTERNAL  PROPERTIES 

Exploratory  coring  at  meter  37  on  cross  section  A,  at  meters  50  and  60  on  cross  section  B,  and 
at  meter  43  on  cross  section  C  (see  Fig.  2)  revealed  no  cavities  of  any  type  and  showed  the  ice  to 
be  sound,  as  indicated  by  augering  resistance.  There  were  no  wide  cracks  or  leads  that  transected 
the  ridge  so  an  actual  profile  could  not  be  examined.  However  Kovacs  (unpublished  data)  has 
examined  cross  sections  of  several  similar  multiyear  ridges  in  the  Beaufort  Sea.  Figure  4  shows  a 
general  view  of  such  a  ridge  which  has  split.  Note  that  the  ice  is  extremely  massive.  Figure  5 
shows  a  detail  of  the  internal  structure  of  another  multiyear  ridge.  The  distribution  and  size  of  the 
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Figure  4.  View  ot  the  sail  ol  a  multiyear  pressure  ridge  located  in  the 
Beautorl  Sea.  Note  the  massive  nature  of  the  ice. 


Figure  5.  Internal  structure  ot  a  multiyear  pressure  ridge  from  the  Beautorl 
Sea.  Note  the  lack  ot  voids  between  the  angular  blocks  composing  the  ridge. 
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T*mptralur«,  C 


Figure  6.  Salinity,  temperature  and  brine  volume  profiles  of  the  ice  at  Meter  37  on  cross 

section  A. 


subangular  ice  blocks  incorporated  into  this  ridge  as  well  as  the  refrozen  water  between  the  blocks 
is  clearly  shown.  The  lack  of  interblock  voids  is  striking.  Note  also  the  dark  plankton  band  within 
t  re  largest  ice  block.  Based  on  this  information,  our  coring  observations  and,  as  mentioned  earlier 

*  Z  *  T  CXperienced  When  shiPs  b^ak  through  multiyear  ridges,  we  suggest  that  ' 

he  ice  blocks  composing  most  multiyear  ridges  are  completely  bonded  together  by  the  freezing  of 
interstitial  water. 

Figures  6,  7  and  8  show  the  salinity,  temperature  and  brine  volume  profiles  ootained  by  coring 
into  the  ridge  Figure  8  also  gives  the  density  profile  of  the  ice  to  a  depth  of  some  3  5m  All 
the  ridge  profiles  show  the  salinity  to  be  virtually  zero  at  the  surface. 

As  shown  in  Figures  6  .nd  7  the  ice  remains  fresh  for  a  depth  of  1  m.  The  near-surface  ice 
was  also  found  to  contain  a  large  number  of  air  bubbles,  suggestive  of  extensive  brine  drainage 
This  increase  in  air  content  is  reflected  in  the  sharp  decrease  in  ice  density  near  the  top  of  the 
ridge  sail  (Fir  8).  b  igures  9  and  10  show  similar  profile  data  obtained  from  the  multiyear  ice 
adjoining  the  ridge.  In  all  of  the  profiles  there  is  a  gradual  increase  in  the  salinity  with  depth. 

At  the  bottom  of  ttie  ridge  the  salinity  averages  roughly  4.5  700  whereas  at  the  bottom  of  the  pre¬ 
sumably  undeformed  multiyear  ice,  the  value  is  very  slightly  higher.  All  the  salinity  profiles  show 
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the  wide  scatter  that  is  particularly  characteristic  of  multiyear  icc  (Schwarzaclier  1959, 
nterstemer  1968)  and  to  a  lesser  extent  of  first-year  ice  (Weeks  and  Lee  1962).  A  point  of  interest 
.5  Ilia,  much  of  the  ,ce  represented  by  the  profiles  shown  Future  10  formed  In  a  mel.  pond  and  Is 

he  sal'  „  T. u'V"  ""****  '»  b°"'  *»  **  surroundins  the  corinp  site  ani 

.!«  saltnuy  profile  ttself.  The  pond  clearly  exhibited  a  pronoaneed  salinity  stratification  prior  to 

The  temperature  profiles  shown  in  Figures  6,  s,  2  and  10  all  show  by  the  curvature  in  their 
uppe(  portions  that  a  gTadual  warming  of  the  ice  cover  is  under  way.  This  trend  is,  however  not 
vet  apparent  in  Figure  7  where  the  thick  snow  cover  formed  in  the  lee  of  the  ridge  has  caused  a 
delav  in  the  arrival  of  the  warming  trend  at  the  ice  surface. 

The  brine  volume  profiles  shown  in  Figures  6,  7  and  8  are  as  expected:  effectively  nil  at  the 

a  d  salmit  *  toward  the  response  to  the  increase  in  the  temperature 

and  salinity  values.  Because  or  tne  time  of  year  (March),  the  ice  temperatures  would  be  expected 

)  »  m  ar  t  icir  annual  minimum  corresponding  to  a  maximum  in  ice  strength  (Weeks  and  Assur  1968). 

ir  we  assume  that  the  overall  ridge  is  in  isostatic  equilibrium  we  can  calculate  a  mean  density 
or  the  ridge  ice  by  determining  the  ratio  of  the  ice  above  sea  level  to  the  ice  be  low  sea lev  The 
now  cov^was  converted  to  an  equivalent  ice  volume  by  assuming  an  average  snow  density  of 
;.nn ,  '  rano  was  1  to  1J'  ^responding  to  a  mean  ice  density  of  0.91  g  c„r3 

ndgefnT teTtosrto  Z  de"Si,y  °f  °'90  K  cm"3*  This  ^ests  that  multiyear’ 

rulgts  may  be  close  to  being  in  overall  isostatic  balance.  Although  a  sample  of  one  can  hardlv  be 
conx.de, od  as  proof  of  this  conjee, me.  i,  ,s  reasonable  ,o  assume  ,„a,  this  wen  m™  be  , he  case 
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Figure  9.  Salitity,  temperature  and  brine  volume  profiles  of  the  ice  at  Meter  89  on  cross  section  B, 
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Figure  10.  Salinity,  temperature  and  brine  volume  profiles  of  the  ice  at  Meter  97  on  cross  section  B. 


Figure  11.  Large  free-floating  ridge  sail  observed  in' the  Beaufort  Sea  near  Camp  200. 
The  highest  point  on  the  sail  is  12.8  m  above  sea  level. 


LARGEST  RIDGE  SAIL 

An  extremely  large  pressure  ridge  was  observed  near  one  of  the  oceanographic  stations  (about 
10  km  from  the  AIDJEX  base  camp).  To  the  best  of  our  knowledge  this  ridge  has  the  highest  free- 
floating  sail  yet  sighted  (and  documented)  in  the  Arctic  Ocean  (12.8  m  above  sea  level).  The  ice 
in  the  ridge  was  first  year  and  the  large  ice  block  located  on  the  top  of  the  ridge  was  2.2  m  thick, 

3  m  wide  and  7  m  long.  The  ridge  had  clearly  been  active  in  the  recent  past  and  a  lead  existed 
along  one  side  of  it  at  the  time  of  our  visit.  The  high  portion  of  the  ridge  was  very  limited  in 
lateral  extent  (=  100  m).  If  we  estimate  the  keel  depth  of  this  ridge  by  using  the  (keel  depth/sail 
height)  ratio  of  4.5  observed  by  Kovacs  (1971)  on  first-year  ridges,  we  obtain  ayalue  of  57  m,  some 
10  m  more  than  the  deepest  previously  observed  keel. 
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CONCLUSIONS 

The  multiyear  rid^e  studied  in  this  report  was  large  and  structurally  massive.  All  the  initial 
voids  between  the  blocks  composing  the  ridge  have  been  subsequently  filled  with  solid  ice.  In 
addition  much  of  the  ice  in  the  ridge  has  a  very  low  salinity.  It  is.  therefore,  reasonable  to  assume 
that  ridges  of  this  type  represent  a  significant  impediment  to  shipping  and  should  be  avoided  if  at 
all  possible.  They  also  represent  a  major  problem  in  the  design  of  potential  offshore  structures 
such  as  shipping  terminals  and  drilling  platforms. 
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by 

W.D.  Hibler,  III,  S.F.  Ackley,  W.F.  Weeks  and  A.  Kovacs 


ABSTRACT 

A  spectral  study  of  the  snow  and  ice  topography  on  a  multiyear  ice  floe  has 
shown  that  the  snow  cover,  although  attenuating  the  roughness  amplitude  of  the 
ice  surface,  does  not  cover  it  completely.  In  general  the  snow  surface  variance 
is  lower  by  a  factor  of  %  to  '/,  than  the  ice  surface  variance.  The  correlation  be¬ 
tween  snow  and  ice  surface  roughness  is  highly  significant  for  long  wavelengths 
(>  8  m),  but  fails  to  be  significant  for  short  wavelengths  (<  4  m).  These  results 
-t?ree  with  what  might  be  expected  intuitively  in  that  long  wavelength  variations 
are  not  masked  appreciably  while  short  wavelength  variations  are  well  hidden. 
Although  the  ice  sheet  as  a  whole  is  in  free-floating,  isostatic  equilibrium, 
pronounced  local  deviations  from  isostatic  equilibrium  are  common.  The  trend 
is  for  icc  drafts  to  deviate  more  than  expected  from  isostasy  for  thin  ice  and  less 
than  expected  for  thick  ice.  Estimates  are  also  made  of  the  n  imber  of  ice  thick¬ 
ness  measurements  required  to  obtain  the  mean  thickness  of  the  multiyear  floe  to 
any  specified  accuracy. 


INTRODUCTION 

Little  work  has  been  done  on  the  correlation  between  the  roughness  of  the  upper  and  lower 
surfaces  of  sea  ice.  At  first  thought  this  might  appear  surprising  inasmuch  as  the  surface  rough¬ 
ness  is  important  in  controlling  the  wind  and  water  drag  exerted  on  the  ice.  However,  there  is 
still  no  indirect  method  of  determining  the  thickness  of  sea  ice.  much  less  a  method  of  doing  this 
remotely  as  a  continuous  profile.  Drilling  is,  of  course,  time  consuming  and  tiring,  but  is  virtually 
the  only  means  of  determining  freeboard  and  draft  at  a  point.  Present  observations  on  ridging  in 
first -year  ice  suggest  that  although  lack  of  local  isostatic  equilibrium  is  common,  a  good  overall 
correlation  between  top  and  bottom  roughness  is  to  be  expected.  Whether  this  is  generally  true  for 
multiyear  ice  has  not  been  pre'  iously  studied. 

The  following  is  a  preliminary  look  at  the  freeboard  to  draft  correlation  problem  in  multiyear 
ice.  In  addition  a  study  of  the  relation  between  the  snow  suiface  and  the  true  ice  surface  of  the 
multiyear  floe  was  completed.  This  is  of  interest  because  lemote  sensing  techniques  using  laser 
profilometry  record  the  elevation  of  the  snow  surface  and  not  the  true  ice  surface. 

The  study  site  was  part  of  a  multiyear  floe  approximately  5  km^  in  area,  located  in  the  Beaufort 
Sea  at  the  site  of  the  1971  AIDJEX  pilot  study  about  550  km  north  of  Tuktoyaktuk,  N.W.T.  An 
aerial  photograph  of  the  study  area  is  shown  in  Figure  1  of  Part  II.  Figure  1  of  Part  III  shows  the 
plan  drawing  of  the  profile  lines  where  1-1  refers  to  the  ridge  shown  in  the  aerial  photograph.  At 
the  time  when  the  photo  was  taken  the  ice  surface  was  effectively  scoured  free  of  the  snow  that  is 
shown  in  Figure  2.  The  structure  of  the  large  multiyear  pressure  ridge  in  the  right  portion  of  the 
figure  was  also  studied  and  is  the  subject  of  another  paper  (Kovacs  et  al.  1972).  The  differences 
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Figure  1.  The  locations  of  the  profile  lines  alone  '  ich  snow  and  ice  elevations  and  random  thick¬ 
ness  measurements  were  taken  (see  Fig.  °art  II  for  airphoto  of  multiyear  floe). 


in  roughness  between  first-  and  multiyear  ice  are  quite  evident  from  the  "freckled”  appearance  of 
the  multiyear  ice  as  shewn  in  the  figure  compared  to  the  smooth  first-year  ice  whose  surface  is 
marked  only  by  snowdrift  features. 

Figure  2  shows  cross-sectional  views  of  the  five  elevation  profiles  that  were  studied.  The 
snow  depth  and  the  ice  elevation  above  sea  level  were  determined  at  1-m  intervals  along  each  pro¬ 
file  by  leveling.  In  addition  ice  thicknesses  were  determined  by  drilling  at  31  locations.  The 
figure  clearly  shows  the  undulating  topography  characteristic  of  multiyear  floes. 


RESULTS 


Snow  surface/ice  surface  corelations 

In  nuKh  remote  sensing  data,  in  particular  in  laser  profiles  and  digitized  roughness  plots  made 
from  aerial  photographs,  the  surface  of  the  snow  is  "seen”  rather  than  the  true  ice  surface.  It  was 
apparent  by  observation  at  the  study  site  that  after  a  storm,  snow  drifting  could  change  the  effective 
surface  sampled  by  remote  sensing.  From  a  visual  examination  of  the  profiles  shown  in  Figure  2 
it  initially  appears  that  the  snow  covers  up  much  of  the  interesting  roughness.  A  more  detailed 
examination  of  the  surface  roughness  indicates  that  this  may  not  be  the  case.  Figure  3  shows  both 
snow  and  ice  spectra  along  all  profiles.  The  spectra  were  calculated  using  a  Hamming  spectral 
window  and  a  maximum  wave  number  of  50.  The  confidence  limits  are  indicated.  Computational 
detail  and  a  discussion  of  the  particular  interpretation  of  the  power  spectra  of  sea  ice  can  be  found 
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Hibler  and  LeSchack  (1972)  respectively. 
For  our  analysis  here,  the  spectra  may 
be  considered  a  plot  of  the  amount  of 
variance  or  roughness  versus  the  frequen¬ 
cy  or  wavelength  of  the  roughness.  The 
plots  shown  in  Figure  3  are  normalized 
to  the  total  variance  so  that  the  area 
under  the  spectral  curve  in  any  given  fre¬ 
quency  band  gives  the  total  variance  or 
mean  square  roughness  in  that  region. 

From  Figure  3  we  see  that  as  ex¬ 
pected  the  spectral  components  for  the 
snow  surface  are  generally  all  smaller 
by  a  factor  of  */,  to  ‘/4  than  the  spectral 
components  of  the  ice  surface.  This 
means  that  the  roughness  has  been 
smoothed  out  and  consequently  wind 
stress  measurements  (Banke  and  Smith 
1971)  might  be  expected  to  change  as  the 
snow  cover  t  hifts.  Taking  the  square 
root  of  the  variance  to  be  some  measure 
of  the  mean  roughness  height  we  see  that 
a  factor  of  2  variation  in  the  mean  rough¬ 
ness  height  can  easily  be  caused  by  shifts 
in  the  snow  cover. 

The  other  point  obvious  from  Figure 
3  is  that  the  general  shapes  of  the  snow 
surface  and  ice  surface  spectra  are  much 
the  same  for  the  same  profile,  with  many 
of  the  same  spectral  peaks  occurring  in 
the  snow  surface  spectrum  as  in  the  ice 

surface  spectrum.  These  peaks  indicate 
hat  the  snow  cover,  although  attenuating 
the  amplitude  of  the  ice  surface,  does  so 
in  such  a  way  that  the  dominant  spectral 
components,  especially  at  low  frequencies 
are  still  observable.  ’ 

-  h  T°  °b,tain  a  more  quantitative  analysis 
of  how  well  the  snow  surface  correlates 


Figure  3.  Snow  elevation  and  ice  eleva¬ 
tion  power  spectra  along  five  prefiles. 

he  spectra  are  normalized  so  that  the 
total  area  under  the  spectral  curves 
equals  the  total  variance  tor  any  given 
profile.  The  y-axis  amplitude  represents 
e  variance  per  frequency  interval.  The 
80, -o  confidence  limits  are  as  indicated. 
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Figure  5.  Unaltered  and  altered  snow  and  ice  profile  data.  The  frequency  response  of  the  band  pass  filters  used  is  given  in  Figure  4.  For  each 
profile,  part  a  represents  the  unfiltered  profiles  and  parts  b,  c  and  d  represent  the  low,  medium  and  high  pass  filtered  results  respectively. 
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the  correlation  coefficient  (using  16  degrees  of  freedom)  is  given  by  the  dotted  line.  The  plot 
clearly  indicates  the  stronger  correlation  for  wavelengths  longer  than  8  m.  In  fact  the  average 
correlation  coefficient  for  the  lowest  frequency  band  is  0.80,  which  is  greater  than  the  average 
total  correlation  coefficient  (for  all  frequencies)  of  0.72. 

The  filtered  profiles  for  each  of  the  three  pass  bands  together  with  ‘.he  unfiltered  results  are 
given  in  Figure  5.  Thise  plots  qualitatively  give  the  same  results  indicated  by  the  correlation 
coefficients,  namely  thit  wavelengths  longer  than  about  8  m  are  well  correlated  whereas  shorter 
wavelengths  have  a  poorer  correlation.  In  summary  the  quantitative  results  suggest  that  as  a  rule 
of  thumb  8  m  is  the  demarcation  wavelength,  with  shower  wavelength  variation  in  the  ice  surface 
being  masked  by  the  snow  cover  while  the  longer  wavelength  components  are  not  strongly  affected. 

Both  before  and  after  these  measurements  were  made,  high  winds  blew  away  much  of  the  snow 
cover.  Consequently,  the  resulting  surface  roughness  spectrum  changed  to  one  more  nearly 
approaching  the  true  ice  surface.  The  spectra  plotted  in  this  paper  give  some  indication  of  the 
limits  of  the  surface  roughness  spectrum  expected  due  to  such  changing  conditions. 

The  spectra  of  the  true  ice  surface  also  indicated  a  dominant  low  frequency  spectral  peak  at  a 
wavelength  of  about  30  m  which  persisted  on  parallel  profiles.  We  feel  this  peak  >s  representative 
of  melt  hummock  spacings  for  the  region  analyzed.  The  low  frequency  peak  did  shift  as  the  profile 
direction  was  changed,  indicating  lineation  (Hibler  and  LeSchack  1972)  in  the  surface  approximately 
parallel  to  the  nearby  pressure  ridge  shown  in  Figure  1.  It  should  be  noted  here  that  the  profiles 
analyzed  do  not  include  the  pressure  ridge,  so  that  we  are  discussing  here  only  the  “flat”  part  of 
a  multiyear  floe. 

Freeboard/ice  thickness  correlations 

To  study  the  correlation  between  surface  roughness  and  thickness  and  the  bottom  topography 
we  calculated  a  correlation  coefficient  of  0.60  for  the  correlation  of  freeboard  to  draft  and  a  co¬ 
efficient  of  0.79  for  the  correlation  of  freeboard  to  thickness.  Both  these  numbers  indicate  a  linear 
correlation  significant  at  better  than  the  1%  level  by  an  analysis  of  variance  test. 

To  determine  the  ice  draft  and/or  the  ice  thickness  from  the  freeboard,  we  can  calculate  the 
regression  line  between  these  variables.  The  resulting  equations  are  d  1.3371  +  2.80  and  t 
2.32/  f  2.80  where  d  and  t  are  the  ice  draft  and  the  thickness  respectively  and  /  is  the  freeboard 
(all  in  meters).  Using  these  relations,  the  root  mean  square  errors  between  the  predicted  and  ob¬ 
served  values  of  draft  and  thickness  were  in  both  cases  0.36  m.  These  numbers  indicate  the 
improvement  that  can  be  made  in  draft  or  thickness  estimates  by  having  inlormation  of  the  local 
freeboard  available.  (The  standard  error  of  the  estimate  of  0.36  m  obtained  by  using  the  regression 
line  is  significantly  less  than  the  standard  errors  of  the  estimate  on  draft  (0.46  m)  and  thickness 
(0.60  m),  see  Table  1). 


h  I.  Summary  of  statistical 

results  (31  holes) 

Mean  ice  draft 

3.24  m 

Ice  draft  standard  deviation 

0.46  m 

Mean  ice  thickness 

3.57  m 

Thickness  standard  deviation 

0.60  m 

Mean  freeboard 

0.33  m 

Freeboard  standard  deviation 

0.20  m 

Mean  snow  depth 

0.16  m 

Snow  deptli  standard  deviation 

0.13  m 

Correlation  coefficients 

Depth  to  freeboard 

0.60  in 

Thickness  ic  freeboard 

0.79  m 
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Figure  6  Ice  freeboard/ice  depth  ratio  (I/d)  versus  ice  thickness  (t).  Thejsostatic  Glance  line 
was  calculated  assuming  ice  and  sea  water  densities  of  0.91  and  1.03  g  cm  .  In  the  effect 
isostatic  balance  the  weight  of  the  snow  was  also  included. 


Both  linear  regression  lines  deviate  significantly  from  the  result  expected  if  the  ice  is  in 
isostatic  balance.  In  particular  the  isostatic  balance  result  would  be  of  the  form  t  =  ad  where  a 
(density  of  sea  water)/(density  of  ice).  Assuming  a  sea  icc  density  of  0.91  g  cm  3 /  e“ed  fr0m 
the  ice  cores)  and  a  sea  water  density  of  1.03  g  cm'3,  we  would  obtain  ratios  for  d  f  and 

7.58  and  8.58  respectively. 

The  trend  of  the  deviations  from  isostatic  balance  is  indicated  very  clearly  in  Figure  6  where 
we  have  plotted  the  freeboard /draft  ratio  vs  thickness.  The  isostatic  balance  line  with  the  snow 
cover  neglected  and  an  effective  isostatic  balance  line  considering  both  the  density  of  the  ice  and 
t,u  dep res  o  of  the  ice  sheet  by  the  snow  cover  are  indicated.  In  calculating  the  effective  line  a 
meansuow  thickness  of  0.16  m  and  a  snow  density  of  0.45  g  cm'3  were  used,  based  on  snow  mea¬ 
surements  at  the  site.  Also  shown  is  the  regression  line  of  f/d  on  t.  The  trend  is  cli  -uly  or 
ice  draft  to  be  greater  than  expected  from  isostatic  balance  for  thin  ice  and  the  inverse  for  thick 
ice,  Therefore,  on  a  point  to  point  basis,  the  assumption  of  isostatic  balance  breaks  down. 

However,  as  might  be  expected  since  the  ice  sheet  is  clearly  free-floating,  if  cross-sectional 
areas  rather  than  individual  points  are  used  for  comparison  with  the  assumption  of  isostatic  balance, 
the  results  are  quite  good.  To  illustrate  this  we  note  that  the  mean  thickness  and  mean  draft  are 
proportional  to  the  total  cross-sectional  area  and  cross-sectional  area  below  the  water  level 
respectively  of  the  ice  floe.  If  a  depression  of  0.07  m  due  to  0.16  m  of  snow  cover  is  included  we 
find  that  the  average  draft  of  3.24  m  and  the  average  ice  thickness  of  3.57  m  imply  a  mean  ice  den¬ 
sity  of  0.915  g  cm'3  for  isostatic  balance.  This  is  in  agreement  with  the  findings  of  Kovacs  et  al. 
(1972)  and  is  a  very  reasonable  value  since  cores  taken  in  this  region  of  the  floe  indicated  densities 
between  0.90  and  0.92  g  cm"3.  The  mean  thickness  we  obtained  also  compares  well  with  the  mean 
value  for  Beaufort  Sea  ''unhummocked”  multiyear  ice  (3.47  m)  obtained  by  Koerner  (1971). 


The  fact  that  individual  points  deviate  significantly  from  isostatic  balance  while  on  the  average 
the  overall  ice  sheet  is  in  isostatic  balance  is  similar  to  the  situation  observed  in  pressure  ridges 
(Fukutomi  and  Kusunoki  1951.  Weeks  et  al.  1972).  The  form  of  the  deviations  suggests  that  a  model 
similar  to  the  Wittinann-Makaroff  pressure  ridge  model  (Wittmann  and  Schule  1966)  may  well  also 
apply  to  the  small  melt  hummocks  that  develop  on  a  multiyear  floe.  In  this  model  the  underside  of 


PART  III.  TOP  AND  BOTTOM  ROUGHNESS  OF  A  MULTIYEAR  ICE  FLOE 


43 


Depth,  meters  Error  (90% confidence  limit)  m 

Figure  7.  Distribution  ol  ice  thickness  Figure  8.  Estimated  error  (90%  conlidence 

measurements.  limits)  ol  the  mean  ice  thickness  as  a  function 

of  the  number  of  holes  used  in  making  the  esti¬ 
mate. 

the  hummock  has  a  greater  horizontal  extent  than  the  sail.  This  would  yield  deviations  of  the  type 
observed,  i.e.  deep  hummocks  are  less  deep  than  would  be  expected  on  the  basis  of  isostasy.  The 
actual  balance  would  be  between  the  volume  of  the  hummock  below  water  and  the  volume  above  water. 

Estimating  the  average  ice  thickness 

One  interesting  application  of  the  thickness  data  is  to  estimate,  using  sampling  theory,  the 
number  of  holes  needed  to  obtain  the  mean  ice  thickness  of  the  multiyear  floe  with  a  specified 
accuracy.  Based  on  the  31  measurements,  the  mean  and  standard  deviation  of  the  ice  thickness  is 
3.57  +  0.60  m.  The  frequency  distribution  of  the  ice  thicknesses  is  shown  in  Figure  7  and  is  to  a 
good  approximation  normal.  The  90%  confidence  limits  on  the  true  mean  are  3.39  and  3.75  m.  If 
we  now  assume  (Griffiths  1967)  that  the  value  of  the  standard  deviation  is  independent  of  the  sample 
size  N ,  a  t  distribution  may  be  used  to  estimate  the  number  of  holes  required  to  estimate  the  ice 
thickness  to  a  given  accuracy.  Figure  8  shows  the  90%  confidence  limits  on  the  mean  calculated 
using  this  procedure  plotted  against  N.  To  obtain  the  mean  thickness  of  the  unridged  portion  of  the 
multiyear  floe  accurate  to  better  than  0.5  m,  with  a  90%  certainty,  only  six  holes  are  necessary.  On 
the  other  hand,  for  an  accuracy  of  0.1  m  with  the  same  certainty,  100  holes  would  be  required. 
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IN  THE  BEAUFORT  SEA 

by 
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ABSTRACT 

“Ice  deformation  in  a  175-km2  area  of  the  Beaufort  Sea  (at  about  74°N,  131°W) 
between  11  and  23  March  1971  was  analyzed  by  means  of  sequential  aerial  photo¬ 
graphs.  The  approximate  time,  direction  and  magnitude  of  deformational  motion 
were  determined  by  combining  observations  of  the  changes  in  features  related  to 
ice  deformation,  such  as  leads,  fractures  and  pressure  ridges,  and  mesoscale 
measurements  of  strains  and  strain  rates  on  the  ice.  During  this  period  a  series 
of  compressions,  extensions  and  shearing  motions  occurred,  producing  numerous 
changes  in  ice  morphology.  The  entire  area  was  apparently  under  considerable 
lateral  constraint  with  no  observable  rotation  or  rounding  or  floes.  No  leads  wider 
than  about  150  m  were  present.  Shearing  motion  along  one  individual  lead  or  up  to 
280  m  was  documented  by  the  photographs.  Relative  motion  along  shear  zones 
was  nearly  always  dextral  or  right  lateral  in  accord  with  the  clockwise  rotation  of 
the  Pacific  Gyre. 


INTRODUCTION 

The  Arctic  Basin  is  dominated  by  a  formidable  cover  of  sea  ice  that  is  in  constant  drifting 
motion.  The  morphology  of  the  ice  is  constantly  changing:  convergence  and  ridging  occur  in  one 
area  while  divergence  with  associated  fracturing  and  floe  separation  occur  in  another  (Kovacs 
1971).  In  order  to  better  understand  the  large-scale  response  of  sea  ice  to  its  environment,  a 
study  of  the  mesoscale  deformation  of  the  pack  ice  in  the  Beaufort  Sea  was  conducted  during 
March  1971  as  part  of  the  1971  Arctic  Ice  Dynamics  Joint  Experiment  (AIDJEX)  pilot  program. 

This  study  consisted  of  measurement  of  strains  and  strain  rates  between  surface  stations  in  a  tri¬ 
angular  network  during  a  two-week  period  from  11  to  23  March  (Hibler  et  al.  1972).  In  addition 
sequential  aerial  photographs  of  the  ice  cover  were  taken  during  this  period  by  NASA  and  U  S 
Navy  planes  at  altitudes  ranging  from  533  to  10,363  m.  This  photography  permitted  observations  of 
the  nature  of  local  ice  deformation. 

This  report  presents  the  results  of  an  analysis  of  some  of  the  air  photos  obtained.  This  study 
is  specifically  relatedto  sea  ice  deformation  in  a  175-km2  area  centered  on  the  strain  triangle. 

Fror  e;uh  set  of  photographs,  overlays  were  prepared  of  features  related  to  deformation  such  as 
leads,  iractures  and  pressure  ridges.  The  overlays  were  then  photographically  reduced  to  a  common 
scale  so  that  the  changes  in  specific  features  could  be  traced  through  the  two-week  period.  When 
used  in  conjunction  with  visual  ground  observations  and  strain  measurements  on  the  ice,  the  approx¬ 
imate  time,  direction  and  magnitude  of  events  associated  with  the  deformation  of  the  sea  ice  could 
be  determined. 
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Figure  l.  Location  map  showing  strain  triangle  and  ice  conditions  on  11  March;  boxed  areas 

are  locations  of  Figures  9-16. 


STUDY  AREA 

The  field  site  for  this  study  was  near  the  Canadian  Polar  Continental  Shelf  Project’s  Camp  200 
in  the  Beaufort  Sea  about  560  km  north  of  Tuktoyaktuk,  N.W.T.,  Canada  (Fig.  1;  see  also  Fig-ire  1 
of  Part  I).  The  camp  was  on  the  edge  of  a  large  multiyear  floe  roughly  16  km  across,  located  at 
73°45'N,  130  15'W.  The  floe  was  surrounded  by  first-year  pack  ice  about  2  m  thick.  This  ice  con¬ 
tained  a  considerable  number  of  pressnre  ridges  and  fractures  as  well  as  some  thinner  zones  of  old 
refrozen  leads.  During  the  two  weeks  of  the  study  the  ice  cover  underwent  considerable  deformation 


METHOD  OF  ANALYSIS 

This  study  is  based  on  aerial  photographs  obtained  during  overflights  by  a  NASA  Convair  990 
and  a  U.S.  Navy  (NAVOCEANO)  NC-121K  on  11,  15,  21  and  23  March  1971  (Table  I).  From  mosaics 
of  these  photographs  overlays  of  1)  leads  and  fractures  and  2)  pressure  ridges  were  prepared  at  a 
common  scale  for  each  of  the  four  days. 

From  these  overlays  changes  in  specific  features  were  traced  through  the  two-week  period  in 
order  to  estimate  the  direction  and  magnitude  of  specific  deformational  events.  For  example, 
extension  or  compression  along  leads  and  fractures  was  determined  by  graphically  restoring  the  ice 
cover  to  its  unfractured  condition  and  measuring  the  approximate  magnitude  and  direction  of  net 
movement  from  the  aerial  photographs.  Like  vise,  shear-type  motions  could  be  easily  recognized 
and  measured  from  observed  offsets  of  pack  ice  features  such  as  old  ridges.  Numerous  new  pressure 
ridges  formed  during  the  period  of  observation  and  in  some  instances  the  amount  of  deformed  ice 
they  contained  could  be  estimated. 

Reference  is  made  here  to  measurements  of  pack  ice  strain  between  the  triangular  network  of 
Stations  a,  fi  and  BC  (base  camp)  established  by  Hibler  et  al.  (1972)  using  distance-measuring 
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Table  I.  Photo  coverage,  1971  AIDJEX  pilot  program. 


Date 

Local 

time* 

Approximate 
altitude,  m 

No.  of 

runs 

Area  covered 

leather  conditions/ 
image  quality 

1 1  March 
(NASA) 

1606-1625 

1,067 

5 

Legs  of  strain  triangle  and  strip 
extending  SE  from  main  camp 

Excellent 

15  March 
(NASA) 

1312-1514 

10.360 

10 

Regional  coverage  N  and  E  of  main 
camp  cut  to  about  35  km 

Fair 

1525-1535 

3.048 

2 

Regional  coverage  NE  of  main  camp 
including  two  strain  stations 

Fair 

1538-1553 

1,067 

3 

Legs  near  but  not  coincident  with 
strain  triangle 

Good 

21  March 
(NAVOCEANO) 

1428-1550 

533 

12 

Strain  triangle  area 

Poor 

23  March 
(NAVOCEANO) 

1355-1417 

529 

3 

Legs  of  strain  triangle 

Excellent 

1448-1650 

1,524 

10 

Regional  coverage 

Good 

‘Local  time  is  7  hours  behind  Greenwich  Mean  Time. 


instruments  (Fig.  1).  Their  data  (Fig.  3)  provided  a  measurement  of  the  net  deformation  (i.e.  strain 
rate  tensor)  for  the  region  while  the  overlays  were  used  to  catalog  and  identify  specific  local  de- 
formational  events  and  history.  The  initial  lengths  of  the  triangle  legs  as  measured  on  12  March 
were:  a  -  fi  8.3  km.  a  -  BC  10.7  km.  and  fi  -  BC  r-  5.6  km. 


ICE  DEFORMATIONS 

11  March 

The  initial  ice  conditions  in  the  study  area  are  shown  in  the  11  March  photographs  (Fig.  1)  and 
in  the  corresponding  overlay  showing  the  lead  and  fracture  system  at  that  time  (Fig.  2).  There  are 
a  large  number  of  existing  ridges.  There  are  also  two  large  leads  (26-35  m  wide)  and  numerous 
smaller  ones  (4-17  m  wide).  Other  discernible  features  are  the  large  refrozen  lead  of  thinner  ice 
near  Station  a  and  the  multiyear  floe  near  BC.  Two  recent  deformational  events  are  clearly  shown 
in  this  photography.  The  first  of  these  is  an  older  localized  NW-SE  shear  of  23  to  29  m  along  three 
fractures  midway  between  a  and  BC.  This  shearing  produced  some  minor  ridges  and  at  the  time  of 
the  photography  the  fractures  had  refrozen  with  a  thin  cover  of  ice  (see  Fig.  lla).  The  second 
event  is  a  major  NE-SW  extension  of  up  to  35  m  that  opened  almost  all  of  the  fractures  in  the  area 
and  was  continuing  at  the  time  of  the  photography  (about  1620).  Strain  measurements  made  between 
1200  and  1820  on  this  same  day  showed  continuing  extension  of  up  to  98  m  along  the  a-fi  line  and 
31  m  along  a-BC  (Fig.  3).  By  late  in  the  day  ground  observers  noted  that  a  large  number  of  leads 
had  opened  in  the  study  area  and  a  large  fracture  had  split  the  multiyear  floe  on  which  the  camp 
was  located,  opening  a  30-m  lead  along  a  pressure  ridge  roughly  200  m  from  their  living  quarters 
(Hibler  et  al.  1972). 


Distonct  Bttwtsn  Stations, 
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^  ♦♦  Compression 

1  ♦  ♦  Extension 


Figure  2.  Lead  and  fracture  system  on  11  March.  Arrows  indicate  type  and  direction 
of  relative  movement;  approximate  magnitude  is  given  in  meters.  Earlier  movement 

indicated  by  parentheses. 


Figure  3.  Measured  extension  and  compression  along  the  triangular  network,  11  to  23  March, 
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Figure  4.  Lead  and  fracture  system  between  12  and  15  March;  note  NE-SW  shear  zones. 


12  March 

Sometime  during  the  night  the  movement  “reversed”  and  compression  occurred  across  the  en¬ 
tire  strain  network  (Fig.  3).  By  morning  most  of  the  leads  had  closed  and  some  new  pressure  ridges 
formed.  A  large  number  of  cracks  and  leads  were  present  in  the  vicinity  of  /S  but  the  station  was 
still  intact  (Hibler  et  al.  1972).  By  1100  the  a-0  and  a-BC  lines  had  shortened  to  their  original 
lengths  and  readings  were  finally  obtained  along  /S-BC,  thereby  "closing”  the  strain  network. 
During  the  afternoon  a-/3  and  a-BC  showed  13  m  and  1  m  compression  respectively  whereas  /S-BC 
had  5  m  extension  (Fig.  3). 

13-14  March 

Although  no  photography  or  strain  measurements  were  obtained  during  this  period,  it  was 
apparent  from  the  strain  measurements,  ground  observations  and  photography  of  15  March  that  addi¬ 
tional  compression  of  up  to  about  35  m  occurred  along  the  a  - /S  and  a-BC  lines  (Fig.  3).  This 
NE-SW  compression  closed  all  of  the  major  fractures  in  the  study  area,  producing  some  ridging  and 
localized  shearing.  In  addition  two  new  E-W  dextral  or  right  lateral*  shear  zones  developed,  one 
across  Station  /3  (7-9  m  shear)  and  the  other  parallel  but  3.5  km  north  of  Station  /S  (up  to  63  m  shear; 
Fig.  4).  This  shearing  motion  was  a  small-scale  replicate  of  a  much  larger  region  J  scale  deforma¬ 
tion  that  was  occurring  across  a  series  of  nearly  parallel  shear  zones,  suggesting  that  the  ice  cover 
was  deforming  as  several  very  large  floes  (Hibler  et  al.  1972). 

15  March 

No  significant  deformation  was  observed  on  this  day.  Along  /3-BC  several  meters  of  com¬ 
pression  occurred,  whereas  the  other  two  began  to  show  a  few  meters  of  extension  (see  Fig.  3). 


♦Movement  of  right  side  down  or  to  the  right  relative  to  the  left  side. 
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o  i  |  l  « 


Figure  5.  Lead  and  fracture  system  be'ween  It 

shear  i 


•*~7~  Shaor 

*  *  Co^ipr»i«lon 

♦  ♦  Extension 

and  20  March;  note  large  extension  to  left  and 
center. 


16-20  Much 

Although  no  photography  or  strain  measurements  were  obtained  during  this  period,  examination 
of  the  15  and  21  March  photographs  and  the  strain  data  of  21  March  indicate  that  extensive  deforma¬ 
tion  occurred  across  the  entire  study  area  (Pig.  5).  Active  E-W  dilatation  apparently  continued, 
opening  two  long  leads,  one  west  of  0-BC  about  142-149  m  and  the  other  near  a  about  41-62  m,  and 
numerous  smaller  ones  (4-22  m).  A  major  dextral  NW-SE  shear  zone  with  56-104  m  movement  de¬ 
veloped  across  the  center  of  the  strain  network.  This  movement  apparently  also  caused  an  8-19  m 
dilatation  along  the  old  shear  zone  running  through  In  addition  a  lead  formed  along  the  pressure 
ridge  at  a,  destroying  the  tent  which  marked  the  station. 

21  March 

By  the  time  of  the  photography  up  to  22  m  of  NE-SW  dilatation  had  occurred  along  three  newly 
refrozen  leads  (Fig.  6).  Along  a-BC  and  a-/3  several  meters  of  continuing  dilatation  was  measured, 
whereas  /9-BC  showed  no  change  (Fig.  3). 

22  March 

No  photography  or  strain  measurements  were  obtained  on  this  day,  but  from  the  photographs  and 
strain  measurements  of  21  and  23  March  it  is  apparent  that  extensive  deformation  occurred  (Fig.  7). 
Major  NW-SE  dilatation  of  up  to  216  m  with  some  shearing  and  new  ridging  occurred  in  the  center  of 
the  study  area  and  east  of  a  (50-168  m).  The  long  lead  west  of  /3-BC  had  74-89  m  of  E-W  com¬ 
pression  but  did  not  close  completely. 

23  March 


During  the  night  a  general  N-S  compression  occurred,  causing  offset  of  47-108  m  in  the  refrozen 
lead  west  of  the  strain  triangle,  69  m  in  the  center  of  the  study  area,  and  36  m  east  of  a  (Fig.  8). 
The  strain  measurements  show  that  very  little  deformation  occuricd  during  the  day  (Fig.  3). 
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RELATIVE  MOVEMENT:  m 
Stoor 

♦  ♦  Competition 

♦  ♦  Extension 


Figure  6.  Lead  and  fracture  system  on  21  March. 


Figure  t .  Lead  and  fracture  system  on  22  March;  note  continued  shearing  in  center  of  strain  tri¬ 
angle. 
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Figure  8.  Lead  and  fracture  system  on  23  March. 


NET  DEFORMATIONAL  CHANGES 

The  deformations  which  occurred  in  this  study  area  over  the  two  weeks  of  observations  are 
photographically  documented  using  comparative  photos  of  11  and  23  March  (see  location  map, 

Fig.  1). 

Station  a 

Extensive  changes  occurred  in  this  area.  In  Figure  9  the  narrow  lead  on  the  right  was  slightly 
compressed,  before  refreezing  (Fig.  9b).  The  new  lead  visible  in  Figure  9  forme  1  on  the  22nd  due 
to  NE-SW  dilatation.  Note  how  it  cuts  across  numerous  pressure  ridges  where  the  ice  is  somewhat 
thicker  and  causes  only  minor  cracking  in  the  relatively  thin  refrozen  lead.  At  Station  a  extensive 
fracturing,  30  m  of  lateral  offset,  and  several  new  pressure  ridges  developed  during  16-20  March 
(Fig.  10).  Fracturing  developed  parallel  to  the  main  ridge  and  helped  produce  a  new  ridge  adjacent 
to  the  original  ridge  (Fig.  10b). 

a-BC  Line 

Along  this  leg  of  the  strain  triangle  leads  closed  from  the  east  and  west  on  13-14  March,  re¬ 
freezing  and  causing  numerous  small  pressure  ridges.  In  Figure  lib  a  20-m-wide  strip  of  ice  was 
compressed  into  a  new  ridge.  Note  how  solid  the  original  fracture  zone  has  become  due  to  refreezing 
and  how  completely  the  original  lead  has  healed  due  to  the  parallel  directions  of  extension  and 
compression.  Figure  12  clearly  shows  the  large  lead  which  developed  on  16-20  March  and  further 
sheared  on  22  March,  resulting  in  about  250  m  of  net  dextral  offset  in  a  NW-SE  direction.  Note  the 
continuing  NW-SE  offset  that  was  occurring  in  the  refrozen  ice  across  the  lead  at  the  time  the  photo 
was  taken  (about  69  m;  Fiv.  12b). 
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a.  11  March.  b.  23  March, 

figure  13.  Ridge  that  formed  near  Station  BC. 


Station  BC 

There  was  relatively  little  deformation  in  this  area.  Figure  13  shows  the  only  major  change, 
a  fracture  which  sheared  about  6  m  NE-SW  during  16*20  March,  creating  a  long,  prominent  ridge. 

Note  its  straightness  and  remarkably  uniform  vertical  profile  as  shown  by  shadows  (Fig.  13b).  This 
new  ridge  looks  very  similar  to  the  much  older  ones  which  surround  it.  The  uplifted  blocks  such  as 
the  one  in  the  upper  right  also  show  the  thickness  of  the  ice  (about  2  m). 

/9-BC  Line 

This  leg  of  the  strain  triangle  showed  only  a  small  amount  of  deformation  during  the  study 
(Fig.  3).  T'te  only  significant  change  was  the  closing  of  the  lead  visible  on  11  March  and  the 
formation  of  a  small  ridge  (Fig.  14b).  This  closure  was  somewhat  unusual  in  that  the  lead  apparently 
refroze  and  then  fractured  again  along  the  southern  edge,  forming  a  small,  irregular  ridge  during  the 
subsequent  compression. 

Station  /3 

A  number  of  changes  occurred  in  this  area.  The  lead  present  on  11  March  (Fig.  15a)  closed  and 
refroze,  forming  a  new  ridge  (Fig.  15b).  The  fracture  visible  on  11  March  also  formed  a  new  ridge 
during  the  same  compressional  event  (13-14  March);  later,  on  16-20  March,  it  refractured  and  the  60- 
m-wide  lead  formed  (Fig.  15b).  Note  the  NE-SW  compression  with  attendant  finger  rafting  in  the 
refrozen  lead  at  the  time  of  photography. 

a-/3  Line 

Along  this  leg  of  the  strain  network  extensive  dextral  shear  occurred,  forming  the  largest  new 
ridge  observed  anywhere  in  the  study  area  (Fig.  16b).  The  shear  motion,  like  that  in  Figure  12, 
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a.  11  March.  />.  23  March. 

Figure  16.  Large  dextral  shear  and  new  ridge  that  formed  along  the  a-  ft  leg. 


started  on  16-20  March  and  further  sheared  on  22  March,  resulting  in  a  net  right-lateral  NW-SE  move¬ 
ment  of  about  280  m.  This  shearing  motion  deformed  a  strip  of  ice  up  to  65  m  wide,  building  a  long, 
nearly  continuous  ridge  about  5-10  m  wide  and  3-4  m  high  (Fig,  16b). 


PRESSURE  RIDGE  DISTRIBUTION 

The  distribution  of  all  pressure  ridges  more  than  about  0.6  m  high  in  the  study  area  as  of 
23  March  is  shown  in  Figure  17.  These  data  clearly  demonstrate  the  ubiquitous  nature  of  the  ridges. 
In  order  to  measure  the  orientation  characteristics  of  these  ridges,  a  representative  study  sector  was 
arbitrarily  chosen  (square  outlined  in  Figure  17).  All  ridges  within  it  were  subdivided  into  segments 
which  could  be  approximated  by  straight  lines.  Next  the  orientation  and  length  of  each  segment 
were  tabulated.  The  mean  segment  length  was  calculated  and  then  used  to  define  a  “standard  ridge" 
(122.1  m).  The  frequency  distribution  of  standard  ridges  as  a  function  of  direction  (Fig.  18)  shows 

a  preferred  orientation  in  the  eastern  and  southern  quadrants  but  a  random  orientation  in  the  northern 
quadrant. 

The  Rayleigh  test  techniques  of  Curray  U956)  and  Krumbein  (1939)  were  applied  to  these  data 
to  define  a  preferred”  orientation  (127°)  and  to  test  if  the  observed  ridges  were  randomly  oriented. 
The  results  yielded  a  probability  of  <  0.00001  that  the  observed  distribution  could  have  been  drawn 
from  a  randomly  oriented  population  (Mock  et  al.  1972).  Similar  results  were  obtained  with  a 
variance  ratio  test  and  a  goodness  of  fit  test  (Mock  et  al.  1972).  From  these  results  and  the  apparent 
absence  of  large  rotational  movements  of  the  ice  "floes,"  we  conclude  that  the  dominant  ridge¬ 
forming  compressions  and  shears  have  come  from  the  east  and  west. 
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Figure  IS.  Frequency  distribution  o(  "stundurd"  ridges  as  a 
function  of  direction. 


SUMMARY  AND  CONCLUSIONS 

During  t he  two  weeks  of  this  study,  at  a  time  when  no  major  cyelonic  passages  occurred,  a 
variety  of  deformation  events  were  observed  in  this  region.  Five  main  types  of  deformation  were 
observed.  1)  fracturing,  2)  extension  along  fractures  and  formation  of  leads  up  to  130  m  wide, 

3)  compression  of  leads  and  crumpling  of  their  thin,  refrozen  ice  cover  into  small  ridges,  4)  lateral 
compression  of  adjacent  ice  sheets  and  creation  of  ridges,  and  3)  shearing  motion  of  up  to  2K0  in 
and  creation  of  ridges.  During  the  study  period  the  entire  area  was  apparently  under  considerable 
lateral  constraint  with  no  significant  rotation  or  rounding  of  floes  such  as  Anderson  (1970)  observed 
in  pack  ice  72  km  northeast  of  Point  Barrow.  Shearing  motion  was  nearly  always  dextral  or  right 
lateral  in  accord  witli  the  clockwise  rotation  of  the  Pacific  Gyre.  At  many  locations  fractures  and 
leads  developed  with  apparent  “disregard”  for  ice  thickness  and  tlnnne-  zones  of  refrozen  leads 
(such  as  in  Figu.es  9  and  10). 

This  study  demonstrates  t hat  a  rather  complete  interpretation  of  the  deformation  history  of  a 
given  sector  of  sea  ice  may  be  made  using  sequential  aerial  photography  in  conjunction  with  surface 
strain  measurements.  For  optimum  results  overlapping  flight  lines  should  be  flown  at  several  dif¬ 
ferent  altitudes.  With  a  13.2-nim  (6-in.)  lens,  the  following  flight  altitudes  are  recommended: 

1.  8000-10,000  m  for  regional  coverage  of  major  lead  systems  and  distribution  of  large  multi¬ 
year  floes. 

2.  1300-2000  m  for  distribution  of  fractures,  leads,  pressure  ridges  and  ice  types. 

3.  600-800  m  for  detailed  delineation  of  deformation  features  for  sequential  studies. 

Even  witli  low-altitude  j  holography  such  as  that  of  21  March  (533  m)  or  23  March  (529  m)  it 
was  impossible  to  detect  any  significant  deformation  in  “passive”  features  away  from  fractures, 
leads  or  relatively  new  pressure  ridges. 
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PART  V:  DATA  ON  MORPHOLOGICAL  AND  PHYSICAL 
CHARACTERISTICS  OF  SEA  ICE  IN  THE  BEAUFORT  SEA 

by 

S.F.  Ackley,  A.  Kovacs,  W.F.  Weeks  and  W.D.  Hibler,  III 


Introduction 

The  data  repined  here  represent  incidental  measurements  made  at  the  main  AIDJEX  Camp  200 
and  at  the  strain  network  camps.  In  general,  the  data  are  not  sufficiently  representative  to  draw 
detailed  conclusions  from  so  the  results  are  listed  only  to  make  them  accessible  and  to  document 
them  for  possible  future  reference. 

Block  size  determinations  in  newly  formed  ridges 

Two  measurements  were  made  of  block  sizes  in  newly  formed  pressure  ridges.  This  parameter 
is  of  some  interest  since  a  recent  model  of  the  ridging  process  in  sea  ice  has  indicated  that  experi- 
mentally*measured  block  sizes  assist  in  the  calculation  of  the  fracture  energy  necessary  to  build 
pressure  ridges  (Parmerter  and  Coon  1973)*. 

The  two  sets  of  measurements  were  made  at  two  locations  near  strain  site  p  (Fig.  1).  The 
ridge  at  site  A  (from  ground  observations)  apparently  formed  by  a  compression  with  the  principal 
compressive  direction  normal  to  the  edge  of  the  thin  ice  in  that  area.  This  ridge  thus  corresponds 
roughly  to  the  case  modeled  by  Parmerter  and  Coon.  If  we  reasonably  assume  that  the  same  direc¬ 
tion  of  motion  caused  the  ridge  at  site  B,  then  this  ridge  formed  by  a  motion  roughly  parallel  to  the 
thin  ice  in  that  area  or  a  shearing  deformation  mode.  The  data  on  block  sizes  for  each  ridge  are 
listed  in  Tat.e  I.  The  block  thickness  was  0.19  m  for  both  ridges. 


Table  I.  Block  sizes. 


Site  A 

Site  B 

Block  width 
(m) 

No.  ol 
blocks 

Block  width 
(m) 

No.  ol 
blocks 

^0.50 

1 

0.90  to  0.95 

3 

0.50  to  0.55 

3 

0.96  to  0.99 

0 

0.56  to  0.60 

5 

1.00  to  1.05 

2 

0.61  to  0.65 

5 

1.05  to  1.1 

0 

0.66  to  0.70 

3 

1.1  to  1.15 

3 

0.71  to  0.75 

5 

1.15  to  1.2 

1 

0.76  to  0.8C 

2 

1.25  to  1.3 

0 

0.81  to  0.85 

1 

1.31  to  1.35 

2 

0.86  to  0.90 

4 

1.36  to  1,4 

1 

0.91  to  0  95 

0 

1.41  to  1.55 

0 

0.96  to  0.99 

0 

1.55  to  1.60 

2 

1  00  to  1.05 

2 

1.20  to  1.25 

2 

•Parmerter,  R.R.  and  M.D.  Coon  (1973)  Mechanical  models  of  ridging  in  the  arctic  sea  ice  cover. 
AIDJEX  Bulletin  No.  19,  Division  of  Marine  Resources,  University  of  Washington. 
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Figure  I.  Location  of  sites  for  block  size  determinations 
approximately  300  m  west  of  Camp  fi.  Data  from  the  two 
sites  are  shown  in  Tables  I  and  II. 

Table  □.  Statistical  averages  of  the  measured  values. 

Site  A  Site  fi 

Mean  block  width  0.73  m  i  .19  m  1.19  m  *  .22  m 

No.  of  measurements  33  14 

Mean  width 'thickness  3.8  6.2 


Salinity  -  temperature  measure’-.o.-its 

In  addition  to  tire  measurements  reported  in  Part  II,  Structure  of  a  Multiyear  Pressure  Ridge, 
two  additional  vertical  profiles  were  also  obtained  (Fig.  2,  3)  as  well  as  some  near-surface  salinities 
in  thick  first-year  ice  and  in  newly  formed  ice  (Table  III)  measured  in  concurrence  with  ground-based 
microwave  ellipsometer  studies  conducted  by  A.  Edgerton  and  D.  Williams  of  Aerojet-General  Corp. 

Two  salinity  values  of  freshly  formed  ice  in  a  “pond”  cut  from  a  new  lead  northwest  of  camp 
for  ellipsometer  studies  are  listed  below.  The  salinities  are  average  values  over  the  thicknesses 
indicated  (sea  water  salinity  32%,,). 

Ice  thickness  Salinity 

3.8  cm  18°/m 

8  cm  207m 

These  sample  were  collected  soon  after  the  ice  formed  and  the  salinity  values  represent  that 
of  newly  formed  ice.  probably  with  very  little  brine  drainage. 
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Temperature,  C  Density,  g/cm3 

a  b 

Figure  2.  Salinity-temperature  profile  from  station  133,  profile  line  C.  (Refer  to  Parts  U  and  III  of 
this  report  for  details  on  the  profile  line  locations.) 


Salinity,  %o 


Figure  3.  Salinity-temperature  profile  from  first-year  ice  located 
just  beyond  the  north  end  of  the  runway  at  the  main  camp.  Ice  thick¬ 
ness  =  1.85  m.  Depth  measured  down  from  top  ice  surf  a Six  cen¬ 
timeters  of  snow  on  top  surface. 

Table  10.  Surface  salinity  and  temperature  measurements. 

Measurements  made  in  same  general  locale  (beyond  north  end  of 
runway  on  first-year  ice)  as  profile  shown  in  Figure  3. 


Temp 

Depth 

Salinity 

<°C) 

(cm) 

<°u 

Site  A 

-19.0 

1 

8 

10.0 

Site  B 

-19.4 

1 

7 

6.4 

